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Abstract 


With the proliferation of modern monitoring and laboratory procedures, physicians 
in intensive care areas may face “information overload”, in dealing with very large, 
complex and ever-changing quantities of clinical data, which often lacks efficient or- 
ganization. ‘This research analyzes the medical knowledge required for formulating 
decision models in the domain of hemodynamics. Based on such analysis, a knowl- 
edge based expert system to track a patient’s hemodynamic state has been developed 
and evaluated in a laboratory setting. 

The initial phase of the work utilizes a cardiovascular simulator to generate “pseudo- 
ICU” waveforms as input to the expert system in order to guide the development of 
the matrix of rules and search strategies. A number of pathological simulations have 
been successfully analyzed by this model-based expert system, including examples of 
hypertension, left ventricular failure, hypovolemia, pulmonary hypertension, etc. We 
conclude that our approach is practical, and provides a mechanism for transforming 
and reducing real-time physiologic data into pathophysiologic hypotheses relevant to 
the management of patients. 
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Chapter 1 


Introduction 


Intensive care units (ICU) provide continuous and often invasive measurements of 
respiratory and hemodynamic status in acutely ill patients. Such monitoring enables 
early detection of changes in the patient’s condition and provides information that 
both directs therapy and assists in evaluating the response to treatment. Needless to 
say, it generates enormous amounts of real-time and off-line data relating to the status 
of these patients. The ICU medical staff must re-assess these patients frequently on 
the basis of data from clinical observations, bedside monitors, mechanical ventilators 
and a wide variety of lab tests. Clinicians usually face a very large, complex and 
ever-changing body of data, and it is often challenging and time consuming for them 
to analyze such massive information loads. As we can imagine, with the development 
of new techniques and medical equipment, the problem of information overloading 
will become more and more serious. The great responsibilities of ICU staff and high 
performance expectations can lead to both emotional and physical fatigue [39]. All 
of these factors may cause errors in patient care. Providing life support in the ICU is 
becoming an increasingly complex task as the volume of monitoring data increases. 
Improvement in the organization and interpretation of the clinical data will have the 
potential of increasing quality and efficiency in the ICU. 

Most intensive care units have the basic capacity to monitor and record heart 
rate and rhythm (usually with a built-in memory and recall capability), venous pres- 


sure, pulmonary arterial pressure, and systemic arterial pressure. In addition, many 
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units have the instruments necessary for measurement of cardiac output. In this 
thesis, I design and implement a knowledge based system which interprets simulated 
ICU hemodynamic data. First, multiple ICU measurable data are gathered from a 
cardiovascular simulator, including central venous pressure, arterial blood pressure, 
pulmonary artery pressure, ventricular pressure, cardiac output, etc. Then I de- 
termine the meaning of these measurements with reference to each pseudo-patient’s 
clinical problem, and suggest appropriate settings for the systemic parameters of the 
cardiovascular simulator to simulate the state of the pseudo-patient. The system is 
intended to help physicians to diagnose and track the disease of patients. 

Thorough analysis of many clinical decision making systems has given us a clearer 
picture of symbolic approaches to medical decision making. Additional requirements 
for decision aids in intensive care areas include programs which can handle data that 
are changing over time in order to provide tracking of the patient’s status during the 
course of the underlying disease or in response to therapeutic intervention. The long 
term objective of this research is to understand how we can manage ICU data in a 
more systematic way, and represent the patients status as a set of parameters in a car- 
diovascular simulator. We seek to develop methods to improve general model-based 
reasoning, and try to find a quantitative or at least a semi-quantitative reasoning 


method. 


1.1 Artificial Intelligence 


The origin of artificial intelligence can be traced back to the earliest days of ma- 
chine computation [16]. In 1843, Ms. A. A. Byron raised the question of whether 
Babbage’s proposed analytical engine, the first programmable computing machine, 
might “think”. Although not called “artificial intelligence” at that time, lots of work 
was done on machine translation beginning in the early 1950s. It was in 1956 that 
artificial intelligence began as a separate aspect of computer science at the famous 
Dartmouth conference [5]. This period was named the prehistory period in [13], be- 


cause few ideas could be experimentally tested. It extended to about 1960 because 
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there was no adequate computers available then. 

Around 1960, the Dawn Age [13] started. People tended to believe that in ten 
years, they could make computers as “smart” as people are. The Dawn Age was 
marked by some limited successes. Patrick Winston named the period from around 
1965 until about 1970 the Dark Period [13], because there was almost no progress. 

Then from about 1970 to 1975, people started to think about how to impress the 
public, industry, etc. That was when some application systems, such as MYCIN, came 
out. It was called Renaissance [13]. This was followed by the Age of Partnerships 
[13], a period when researchers started to pay attention to the research fruits coming 
out from other domains and to form some important liaisons. 

Then from 1980 till about 1990, it was Age of the Entrepreneur [13]. During that 
period, lots of start-up companies focused on designing possible intelligent systems. 
After 1990, with the failure of almost all of those intelligent-wanna-be companies, I 
would like to call our present age the Age of Adjustment. The name comes from my 
analogy to the stock market. We cannot expect that a tree can grow to the sky. The 
potential and the expectation have been oversold by the media and the public. Now 
it is time to come back to re-adjust our positions and expectations for both long term 
and short term. 

There has been tremendous progress in some of the subfields of artificial intelli- 
gence during the past several decades. Some successes are very domain specific, such 
as: XCON (or R1) developed by Digital Equipment Corporation and Carnegie Mellon 
University for doing computer configuration; PUFF developed by Stanford University 
for analyzing certain lung problems. Other successes are less domain specific, such as 
INTELLECT, a powerful interface between decision makers and databases. 

But the problem is just what Gerald Sussman et al. said [3], “The field overall ... 
is moving increasing attention to near-term applications, retarding progress toward 
comprehensive theories and deep scientific understanding, and ultimately, retarding 
progress toward developing the science needed for higher-impact applications.” Part 
of the reason for people to focus on near-term applications is the growing anxiety 


and lack of trust in AI research among the public. As a result, researchers are very 
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eager to show some progress to the public to make them feel more comfortable and 
confident about artificial intelligence. If 1 may make another similar analogy, we can 
say that the share holders (public) are too desperate about artificial intelligence. But 
we can not blame the general public too much, because I think part of the reason 
for people to have such high expectations is that some researchers have made and 
continue to make rash predictions of rapid progress in the near future. 

An interesting criticism comes from some neuroscientists who claim that computer 
scientists want to make a machine, but they don’t even know how the machine works. 
In their opinion, most artificial intelligence researchers should either find new jobs 
or go home to sleep and wait until they figure out the mechanisms of the brain. An 
easy counter-argument is that we did not have to understand aerodynamics to make 
a plane. 

Other criticisms come from some philosophers who claim that one cannot make 
another equally intelligent system with oneself, just like one cannot lift oneself, no 
matter how strong s/he is. 

Personally, I am optimistic about artificial intelligence in the long term, but pes- 
simistic in the short term. Some people even express the concern that. artificial 
intelligence is the biggest joke in the 21st century. The history of the development of 
science tells us that AI is very unlikely to be a joke. Now it is time for us to focus 
on solving some fundamental questions of artificial intelligence. Shimon Ullman et 
al. said [3], “Traditional approaches to artificial intelligence have been largely orga- 
nized around unconstrained algorithms or simplistic biological models. Both lack of 
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constraints and incorrect constraints have imposed show-stopping penalties.” Many 
biologists believe that the neural network in computer scientists’ dictionary should be 
erased or changed, because that is an over-simplified model for neurons. Then should 
we all go back to the bench and look for a more complicated biological model? It 
is easy to measure the detailed activity of a single neuron, but measuring the sys- 
temic activity of a neural region is very difficult because we lack the needed powerful 


computer hardware and bio-techniques. One thing I can be sure of is that we need 


concrete, not superficial, cooperation between scientists in different fields. We need 
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some team work. It is very naive to try to form some national team to cope with this 
problem, because after all this is not the same case as fighting some epidemic disease 
like AIDS. But forming some regional or institutional cooperative teams is certainly 
feasible. Gerald Sussman, Patrick Winston, Shimon Ullman Kenneth Yip, and other 
colleagues have recently formed a team to try to redirect our research toward a pro- 
gram focused on understanding natural intelligence and building integrated intelligent 
systems. To conclude my thoughts on prospects for artificial intelligence, I will say: 


“The light is on the way, but far away.” 


1.2 Knowledge-Based System 


There are three major research issues concerning the knowledge based paradigm [6]: 


e Knowledge Representation: In which way should knowledge be represented as 


symbolic data structures for computer use? 
e Knowledge Utilization: What designs are available for the inference procedure? 


e Knowledge Acquisition: How can we develop a systematic way for computers 


to access the knowledge? 


The earliest knowledge based systems date back to the mid 1960’s [35]. One of the 
conspicuous systems was MACSYMA, developed by the Mathlab Group at Project 
MAC at M.LT. MACSYMA can be used to provide a very broad range of symbolic 
computing capabilities, such as integration, power series and support for manipu- 
lating matrices [35]. Another major early expert system is DENDRAL, developed 
by Stanford University and used to determine the three dimensional structure of a 
chemical compound from its chemical formula and some other information. 

The success of those early expert systems boosted the development of other expert 
systems. In [10], Feigenbaum reported that there were about 1500 expert system 
applications in actual use and a few thousand more in field testing. Despite the 
broad range of applications and different domains, these systems focus on several 


central tasks [35]: 
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e Diagnosis: Focus on how to get a reasonable explanation for some abnormal 


values of observables. Examples include medical diagnosis, etc. 


e Selection: Focus on how to make intelligent selection. Examples include insur- 


ance, credit authorization, etc. 


e Configuration and planning: Focus on how to complete a partial design or adjust 


a complete design to meet some constraints. Examples include XCON. 


Parallel to the development of artificial intelligence, knowledge-based systems are 


also facing some fundamental problems [35]: 


e Fragility of Encoded Knowledge: How to incorporate common sense knowl- 
edge into the system? How to avoid common-sense-violating behavior? One of 
the biggest issues here is that human beings have too much common sense to 


explicitly encode in computers. 


e Complexity in the Real World: Many systems work well for straightforward 


problems but fail on the more complicated ones. 


e Sources of Better Models: In many domains, we lack the ability to capture or 
design the depth and breadth of the domain knowledge needed to allow the 


computer to produce expert behavior. 


Reviewing these problems, we may feel that we actually know too little about 
ourselves, about our learning process, about our knowledge storage forms, etc. Do 
we believe that the storage form of the knowledge in our brain is kind of synapse 
connection, chemical compound diffusion or something else? One may argue that 
probably the reason why the computer is having problems with common sense and 
complexity is that the computer was made too fast to use. Human beings develop 
common sense after living in society for about 20 years. Imagine a scenario, in which 
we put a computer with a new-born kid together, and each time when we teach this 
kid something, we will try to send the same knowledge to the computer. Would this 


be enough to develop an intelligent machine? But the major problem here is that 
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most of the time, we will have trouble teaching computer to learn something new and 
abstract. Typically for the common sense knowledge, we do not know how to extract 


the features of the common sense knowledge. 


1.3 Artificial Intelligence in Medicine 


It has long been recognized that computers could be put to beneficial use in the 
medical field. As early as the mid 1950’s, physicians and computer scientists had 
recognized that computers could be used to assist in making clinical decisions and 
potentially to implement automatic decision making systems. There are several gen- 


eral approaches [13]: 


e Flowcharts: Encode a serial description of actions which good physicians will 


take. 


e Large Clinical Data Base: Search among previously studied cases and match 


the current case. 


e Probability Theory: Based on the observable data, derive the probabilities of 


several explanations. 


But starting in the 1970’s, scientists involved in computer-based reasoning began 
to recognize the potential benefits of applying symbolic reasoning techniques in clini- 
cal domains [33]. The first medical reasoning program, known as the MYCIN System 
[32], adopted symbolic processing techniques largely in response to a conviction that 
computer-based consultation systems, in order to be accepted by physicians, should 
be able to explain how and why a particular conclusion had been derived [32]. Sub- 
sequently, a series of additional medical application programs have been created. For 
example, the PIP [28] (Present Illness Program) is designed to simulate human clin- 
ical cognition and the INTERNIST [27] diagnostic system of Pople and Myers is a 
computerized diagnostic program which emphasizes a very broad coverage of clinical 


diagnostic situations. These systems were impressive, but their use revealed many 
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more problems, such as inadequate criteria for deciding when a diagnosis is complete 
[36], and the need for better representation of anatomical and physiological medical 
knowledge [17]. 

Although several expert systems were developed for medical problem solving dur- 
ing the past couple of decades, AI has failed to make a significant impact on health 
care delivery [37]. First of all, most initial efforts in artificial intelligence in medicine 
(AIM) were focused on developing automated diagnosis systems, which make medical 
professionals nervous and unwilling to cooperate [37]. Other reasons are the com- 
plexity of medicine itself, including the complexity and inherent variability of human 
anatomy and physiology, and our lack of understanding of the medical decision pro- 
cess on a cognitive level. With the general dim outlook of artificial intelligence, AIM 
is also facing its greatest challenge. We must make some applicable systems to please 
society. Now AIM is focusing to a greater extent on knowledge management and 
intelligent advice-giving instead of the ambitious task of automated diagnosis. It is 
much easier to generate some advice than to replace physicians, both scientifically 
and socially. 

In the past several years, there has been a rise in AIM activity. Some of the 
reasons for the increased activity are the increase in computer power and a rise in 
physicians’ interests in computer technology. Due to the increasing demands of infor- 
mation processing that continue to challenge physicians, they have come to demand 
computer technology for patient care purposes. Therefore, the development of meth- 
ods to manage endless streams of data and to assist in making decisions using these 
vast sources of information are among the top priorities of AIM researchers [32]. On 
the other hand, significant progress has been made in model-based reasoning, mainly 
due to the development of theories and methods for qualitative modeling of physical 
devices and systems. Several model-based computer programs have been devised in 
medical AI, such as AI/MM [20] (dealing with a mathematical model of renal physi- 
ology), LOCALIZE [11] (dealing with the localization of lesions within the peripheral 
nervous system), The Heart Failure Program [21] (reasoning about the function of 


the cardiovascular system), etc. 
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The current problem of AIM is that with our continued lack of fundamental under- 
standing of the reasoning process, how far can AIM go? If we know that intelligence 
can be represented or replaced by pure calculation, then we believe, theoretically, 
that computers can replace or simulate the thinking process of human beings. Even 
so, applying artificial intelligence techniques to the medical system would bring some 


extra concerns which do not exist in other domains: 


e Cost of the Computation: Obviously, we must set a time limit for running 
computers when we need real-time responses. For example, we can not wait a 


day for the computer’s results in ICU areas. 


e Legitimization: How to let those intelligent pseudo-doctors go through the le- 


gitimization process? If the wrong diagnosis happens, who will be responsible? 


e Social Acceptance: Patients may just refuse to be seen by those intelligent 
pseudo-doctors, because they feel more comfortable in dealing with physicians 


or nurses personally. 


1.4 Background and Motivations 


Twenty years ago, ICU care focused solely on the assessment of vital signs and lab- 
oratory values. Today, routine intensive care involves a great variety of monitoring 
devices, infusion pumps, drug-administration systems, sophisticated ventilators, etc. 
The responsibilities of ICU staffs are complex [31]; they usually assess more than 50 
measurements, laboratory values, and physical findings, etc. Also, most of the current 
medical monitoring systems include built-in alarms, and permit users to change the 
alarm limits. However, these automated alarms operate independently, and are so 
unaware of the clinical context that they have a high rate of false alarms. ICU staffs 
may delay their response to alarms that are frequently in error, or even ignore them 
[24]. 

The large amount of monitored ICU data may create information overload, which 


leads to errors and mishaps in ICU care. There are occasional tragedies, most of which 
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are due to human error, reported in ICU care. In [1], Abramson estimated that about 
90 of 150 adverse occurrences in a surgical ICU were caused by human error. We may 
also find the similar results regarding human errors which happened in ICU care 
[29]. So we may expect that with the development of medical instruments, more and 
more clinical data are going to be presented to ICU staffs. So information overload 
will become an increasingly serious problem unless better methods are designed to 
organize and present the data to the clinicians. Relative to the rapid development of 
medical instruments and tests, the slow progress of information process and synthesis 
may impede the overall advancement of health care. Fortunately, computer based 
access to ICU data has made it possible for researchers to analyze large databases 
and to create statistical prediction models which can be used to qualitatively evaluate 
and predict the status of patients, such as APACHE (Acute Physiology, Age, Chronic 
Health Evaluation), which predicts short-term survival in the ICU. 

So far, most of the intelligent alarm systems have failed to interpret the meaning of 
changes of multiple physiologic variables for the overall clinical context. Those alarms 
usually only highlight important events. To perform an overall clinical context-based 
interpretation will require a model for ICU patients, so that we can compare the 
current status of patients with respect to their previous status in an overall clinical 


context, instead of only focusing on one or two physiological signals. 


1.5 Related Work 


In the early 1980’s, Larry Fagan developed an expert system called VM (Ventilator 
Manager). VM is designed to interpret on-line quantitative data in the ICU [9]. It 
was able to produce interpretation of physiological measurements over time. Since 
then, a number of approaches have been developed to interpret ICU data [12]. In 
intensive care areas, automated methods that monitor the hemodynamic status of 
patients must reason about and respond to changes in large amounts of clinical data, 
such as blood pressure, heart rate, cardiac output, etc. Many researchers have devel- 


oped techniques to interpret ICU data [12]. Generally, we can divide these techniques 
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into two categories: numeric and symbolic methods. The numeric methods are rea- 
soning processes for providing quantitative analysis; and symbolic methods deal with 
qualitative analysis. In intensive care areas, the clinical context includes data that are 
both numeric (such as that the arterial blood pressure is 120/80 mmHg) and symbolic 
(such as the interpretation of an echocardiogram). Other programs included qualita- 
tive information in the interpretation of patient data by applying rule-based symbolic 
models. The GUARDIAN project, for example, implements an adaptive intelligent 
agent to reason about monitored ICU data within the real-time constraints of ICU 
care. Mathematical models are also powerful tools for simulating the quantitative 
time-dependent behavior of complex, dynamic physiologic systems. Mathematical 
models that include physiologic concepts and feedback mechanisms allow researchers 
to study interactions among physiologic systems. 

We believe that qualitative and semi-quantitative analyses for patients in the 
ICU are very important when we try to estimate the parameters that describe a 
patient’s physiologic state. We developed a method for finding patient-specific model 
parameters of detailed mathematical models. The patient-specific model parameters 
will be presented to the physicians to illuminate the fundamental pathology and to 


suggest patient-specific treatment recommendations. 


1.6 Overview of Thesis 


This introductory chapter has briefly discussed the development and the current sta- 
tus of artificial intelligence, knowledge-based systems and the development of the spe- 
cific domain — artificial intelligence in medicine. We talked about the background 
and the motivations of this project, and also gave some background about related 
work. 


The remainder of this thesis is organized as follows: 


e Chapter 2 describes the medical domain associated with this work, that is ICU 
patient care. Specifically, it describes some cardiovascular patients. An example 


case is also included for illustration. 
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Chapter 3 describes the tool, a cardiovascular simulator, used in the research. 
It gives an overall picture about the structure of the hemodynamic model and 
baroreflex model, which are used in the simulator. Some example cases to 


simulate cardiovascular patients are also included. 


Chapter 4 details the analysis of the reasoning model. First, we talk about the 
homotopy method, which is the method we used when we started the research. 


Then the method used in our current system is elaborated. 


Chapter 5 describes the evaluation of the system. What is our ideal system? 
How do we evaluate its performance? Finally, we lead to the discussion of the 


question — What is acceptable performance? 


Chapter 6 lists the accomplishments, limitations, practical considerations and 


future directions of this work. 


Ze 


Chapter 2 


Domain Knowledge 


An intensive care unit is defined by its ability to provide the enviroment, facilities, 
and personnel for the care of severely ill patients. The important features required of 


such a unit are listed as follows [39]: 


1. High nurse-to-patient ratio 

2. Ready accessibility of physicians 

3. Ability to provide invasive hemodynamic monitoring 
4, Availability of respiratory support techniques 


5. Ability to provide supervised intravenous infusions of pharmacologic agents 


In severely ill patients, the interrelationship between different physiological sys- 
tems must be noticed. ICU staff should not focus on one component of the illness, 
even if that component is predominant and may yield a therapeutic approach. Usu- 
ally, management of patients in the ICU presents lots of dilemmas. For example, 
treatment directed to reduce the pulmonary capillary wedge pressure in patients with 
adult respiratory distress syndrome may adversely affect renal and central nervous 
system perfusion. On the other hand, increasing the pulmonary capillary wedge pres- 


sure to optimize cardiac output in a patient with ischemic heart disease may result in 
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non-cardiogenic pulmonary edema. In most cases, ICU staff must synthesize an over- 
all diagnostic and treatment strategy that considers the views of various specialists 


with a narrower focus and then make some compromise. 


2.1 ICU Patient Care 


Most ICU areas have a very general orientation, treating almost all types of severely ill 
patients. Because of the sickness of ICU patients, ICU units require clear delineation 
of administrative and medical guidelines of authority and responsibility. Therefore, 
a general guideline for admission and discharge of patients is required. The exis- 
tence of such policies reduces the apparent ambiguity often inherent in the difficult 
environment of a ICU and enables prompt decision making by ICU staff [39]. 

An ICU usually can provide continuous and often invasive measurements of respi- 
ratory and hemodynamic status for severely ill patients. The complexity of monitoring 
systems varies considerably, it can range from simple electrocardiographic monitoring 
with only a real-time screen display to automated “closed loop” systems to regulate 
intravenous infusions of fluids and drugs. 

The complications of intensive care usually are difficult to detect from the illness 
[39]. For example, cardiac arrhythmias and gastrointestinal hemorrhage are very 
common in ICU patients, but they may or may not necessarily relate to the care in 
the ICU. However, some complications are straightforward and clearly related to an 
intervention taking place in an ICU. For example, ventricular tachycardia (VT) can 
occur during the passage of a Swan-Ganz catheter through the right ventricle. 

In addition to patient-related complications, the ICU care environment can also 
stress ICU staff psychologically. There have been many reports that the negative 
psychologic effects of working under great high pressure conditions can cause ICU 


staff “burn-out”. 
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2.2 Respiratory Monitoring 


Usually, as a minimum, respiratory monitoring should include measurement of res- 


piratory rate and periodic measurement of pressure of oxygen (P,,.), pressure of 


005 
carbon dioxide (Pago, ), and arterial blood pH. Respiratory rate can be measured 
and recorded automatically in non-intubated patients using impedance devices. Both 
respiratory rate and tidal volume delivered by the ventilator (Vr) can be monitored 
in intubated, spontaneously breathing patients using a spirometer and appropriate 
alarms [39]. 

For the mechanically ventilated patients, it is essential to monitor respiratory 
rate, exhaled Vr, and airway pressure. We may also use other monitoring techniques, 
such as breath by breath measurements of respiratory system compliance, volume- 
pressure, volume-flow relationships, and measurements of Fig, and exhaled carbon 
dioxide and oxygen. These measurements may provide some useful information such 
as early indications of changes in Py,9,, but their usefulness remains uncertain. 

The transcutaneous Po, and Peo, that are indirect reflections of Fads and Lees 
can be measured continuously using heat skin electrodes. The measurement of Po, 
remains uncertain for adults. Pulse oximeters are noninvasive monitoring devices 
that measure and record oxyhemoglobin saturation. This measurement is useful in 
some clinical situations, such as evaluating oxygenation during sleep and during bron- 


choscopy procedures. Indwelling catheter electrodes for measurement of P,,, and 


PB 


co, and arterial pH have been used but have some technical limitations [39]. 


There are four major complications related to respiratory support in ICU areas 
[39]: 
1. Oxygen Therapy is administered by external devices. It may cause the fol- 


lowing side effects: 


e discomfort related to the device or to administration of dry gas, which can 


be managed by changing the device or improving the humidification 


e fires and hypoventilation because of uncontrolled administration of oxygen 


in patients 


30 


2. Artificial Airways are associated with a number of complications. The im- 
mediate and long term injury to nose, hypopharynx, larynx or trachea, can be 


caused by an endotracheal tube. It may also cause pulmonary infection. 


3. Tracheostomy avoids the problems of an endotracheal tube because it bypasses 
the upper airway. However, tracheostomy has more disadvantages than advan- 
tages. Problems can include hemorrhage, mediastinal and tracheal stenosis. It 


may also cause pulmonary infection. 


4. Mechanical Ventilation can have overventilation and underventilation prob- 
lems. Oxygen toxicity usually occurs in patients being mechanically ventilated 


with gas mixtures containing high concentrations of oxygen [39]. 


2.3. Hemodynamic Monitoring 


Since our project focuses on the interpretation of hemodynamic data, a detailed 
description of hemodynamic monitoring in the ICU is worthwhile. Most ICUs have 
the basic capacity to monitor and record heart rate and rhythm, venous pressure, 
pulmonary arterial pressure, and systemic arterial pressure. Many units also have 
the instruments to measure cardiac output. 


Clinically, the following monitoring is usually done in ICU areas [39]: 


1. Electrocardiagraphic Monitoring is of great value in patients with specific 
cardiac disorders such as acute myocardial infarction or cardiac arrhythmia. 
Also, because it is easy and relatively inexpensive, ECG monitoring has been 
widely used in ICU areas. The automated detection of heart rate and rhythm is 
generally an integral part of the monitoring systems. Abnormalities in heart rate 
or rhythm may signal some worsening of respiration, electrolyte abnormalities, 


and other noncardiac problems. 


2. Systemic Arterial Pressure Monitoring is obviously of great value. Changes 


in blood pressure can be detected immediately and therefore can enable beat by 
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beat assessment of the effects associated with the changes in ventilatory pattern. 
Note that beat to beat variations in systemic arterial blood pressure may not 
warrant specific intervention. The usual technique to monitor systemic arterial 
blood pressure is to insert percutaneously a Teflon catheter into an accessible 
artery. The catheter is connected through a stopcock to a rigid connecting tube 


that in turn is attached to a transducer. 


3. Central Venous Pressure Monitoring is useful in assessing the course of 
right ventricular failure, right ventricular infarction and tricuspid regurgitation. 
The usual technique to monitor central venous pressure is to insert percuta- 
neously a catheter into either the subclavian or the external or internal jugular 


vein. 


4. Pulmonary Arterial Pressure can be measured by a Swan-Ganz catheter. 
This catheter has a balloon just proximal to the tip and a separate lumen 
for inflating the balloon. The Swan-Ganz pulmonary catheter can provide mea- 
surement of pulmonary arterial blood pressure (Pp,), pulmonary arterial wedge 
pressure (Ppaw), central venous pressure (CVP) and cardiac output. However, 
there are some problems with the measurement of Pp,. Because of the oscilla- 
tions in pleural pressure on the measured intravascular pressure, the Pp, signal 


must be carefully analyzed in order to obtain an accurate value. 


5. Cardiac Output can be obtained using the thermistor equipped Swan-Ganz 
catheter. In addition, the instrumentation for performing thermal dilution pro- 
cess and a computer for processing dilution measurements and calculating car- 


diac output are required. 


Based on the above measurements, ICU staff then use Table 2.1 to guide an 
understanding of a patient’s situation. 

There are three principal types of problems associated with hemodynamic mon- 
itoring: local complications related to vascular access, passage and final positioning 


of the catheter, and inappropriate decision making based on inaccurate data or in- 
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Table 2.1: Patterns of Hemodynamic Abnormalities in Severely Ill Patients, from 
[39].Ps4: Mean Systemic Arterial Pressure; Pr4: Mean Right Atrial Pressure; Pp: 
Mean Pulmonary Arterial Pressure; Ppyw: Mean Arterial Wedge Pressure; C.O.: 
Cardiac Output; PVR: Pulmonary Vascular Resistance; SVR: Systemic Vascular 
Resistance 


complete data or misinterpretation of information from the monitoring device. The 


usual complications related to hemodynamic monitoring are listed as follows: 


e Systemic Arterial Pressure: The most frequent complication related to SAP 
monitoring is the formation of a hematoma at the site of the arterial puncture. 


The hematoma may result in the peripheral nerve damage. 


e Central Venous Pressure: The complications of CVP monitoring include air 


embolism, infection and venous thrombosis. 


e Pulmonary Arterial Pressure: All of the complications associated with CVP 
monitoring can happen here. In addition, PAP monitoring can cause a distur- 
bance of cardiac rhythm or conduction or both. Premature ventricular contrac- 


tions occur very commonly when the catheter passes through the right ventricle. 


2.4 An Example Case 


An example case in the domain of ICU patients with severe cardiovascular problems 
is shown below [7]: 
The patient is a 37 year old woman with a positive history for HIV and intravenous 


(IV) drug abuse. She was brought into an ICU with a fever of 104F, extreme lethargy, 
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scant urine output, and a tender abdomen, especially in the right upper quadrant. 
She was noted to be tachycardic, hypotensive, and tachypneic. After being accepted 
to the ICU, her systemic arterial pressure, central venous pressure, right ventricular 
pressure, pulmonary arterial pressure, and pulmonary arterial wedge pressure were 
continuously recorded for assessment. 

In the following chapters, we will examine the process for automated diagnosis and 
propose an iterative method for setting the systemic parameters of a cardiovascular 


simulator to simulate the status of ICU patients. 
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Chapter 3 


Cardiovascular Simulator 


Automated methods that provide treatment advice require a model of physiology to 
assist them in the interpretation of ICU patient data. A lumped-parameter cardio- 
vascular model is used in our expert system to represent the patient’s hemodynamic 
state. 

The cardiovascular simulator (CVsim) is a dynamic computer simulator of hu- 
man cardiovascular hemodynamics [7], originally intended for students of physiology 
and medicine in Harvard-MIT Division of Health Sciences & Technology (HST). It 
is implemented on SUN workstations running the X Window System, allowing stu- 
dents to perform a variety of measurements in a real-time simulated cardiovascular 
system, not all of which would be possible in an animal laboratory. The simulation 
is based on the lumped-parameter mathematical model taught in the Cardiovascular 
Pathophysiology course (HST090) in HST. 

The model is shown in Figure 3-1. It includes four major components: left heart, 
systemic circulation, right heart, and pulmonary circulation. Each side of the heart is 
modeled by a variable capacitor (representing the pumping action of both atrium and 
ventricle), two diodes representing the AV and arterial valves, and outflow resistances. 
There are 23 parameters for defining the status of the simulator; their normal values 


are shown in Table 3.1. 
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Figure 3-1: Circuit Diagram Equivalent of Lumped Parameter Model 
3.1 Background 


Several cardiovascular system models have been put forward since 1950’s [23]. In 
1959, the first model of the cardiovascular system was formulated by Grodins [14]. In 
that model, he defined the steady state operation of the ventricle, then combined two 
ventricles with a pulmonary circuit branch to simulate the effects of Starling heart- 
lung relationship. This work was considered as the starting point for integrative 
physiology [7]. Then in early 1960’s, a six node electrical circuit model was proposed 
by Defares and his colleagues [8]. In that model, each node consists of a capacitor 
to ground and is connected to two adjacent nodes by impedances. The varied left 
and right ventricular capacitances were used to simulate the time varying elastance 


property of contraction. Early in 1970’s, Dr. Arthur Guyton developed refined com- 
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Table 3.1: Normal Values of the Input Parameters in CVsim 


puter models of the cardiovascular system, and a graphical approach for analyzing 
the regulation of the cardiac output based on the model [15]. 

Historically, there were several hemodynamic models designed mainly for teach- 
ing. In [2], Beeuwkes first used the analog computer in the hemodynamics teaching 
laboratory at Harvard Medical School. In 1978, Katz and his colleagues introduced a 
real time digital computer model for hemodynamics [18]. In his model, he simulated 
a Windkessel circulation, which is an over-simplified model in hemodynamic analysis. 


In 1982, Campbell et al. [4] developed a five component model for hemodynamics. In 
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that model, he incorporated linear vascular elastances , resistances and time-varying 
ventricular capacitances proposed in [34]. 

In 1984, as part of the Athena project at M.I.T., Prof. Roger Mark proposed to 
design a cardiovascular model as a teaching tool for Quantitative Physiology and Car- 
diovascular Pathophysiology courses to provide students a quantitative understanding 
of the hemodynamics system. Based on Robert Sah’s model shown in Figure 3-1, 


Timothy Davis implemented it the X window system in 1990 as his Master’s project. 


3.2 Hemodynamic Model 


In analyzing the function of the cardiovascular system, it is convenient to represent 
the hydraulic properties of the various elements in terms of electrical circuit analogs. 


Table 3.2 lists the hydraulic variables and their electrical analogs. 


Hydraulic Variable | Electrical Variable 


Pressure (P) Voltage (V) 
Flow (J) Current (i) 


Volume (Q) Charge (q) 


Capacitance C=22 | Capacitance C=! 


Table 3.2: Analog between Hydraulic System and Electrical System 


First, let us focus on the peripheral circulation. We may characterize the large 
arterial vessels as a single lumped capacitance element, C,. The microvasculature, 
which includes arterioles, capillaries and venules, can be characterized as a single 
resistance R,. The main venous branches can be represented by a lumped venous 
capacitance C’,. In addition, we use a resistance R, to simulate venous return. Then 
we get a simplified electrical circuit model for peripheral vasculature as shown in 
Figure 3-2. Actually, the venous resistance is distributed throughout the whole 
venous system. So a more accurate model of the venous system would be a distributed 


system for the capacitances and resistances. 
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Figure 3-2: Electrical Circuit Model of Peripheral Vasculature 


Now let us consider a model for the heart-lung unit. Then later on, we can 
integrate this unit with the peripheral system to form the model of the whole car- 
diovascular system. In [34], Sagawa et al. stated that if we examined a series of 
ventricular pressure-volume loops for a range of different preloads and afterloads, we 
would find that the end systolic points all lie on a straight line. The diastolic pres- 
sure volume relation may also be described by a straight line. Now we may use a 
variable capacitor to simulate the behavior of left ventricle. Considering the blood 
flow direction in the heart and the function of the cardiac values, we may add two 
diodes in the circuit. In addition, an output resistance is added to the left ventricle 
to simulate the impedance of the blood passing through the aortic valve. Then we 


get the electrical circuit model of the left heart as shown in Figure 3-3. 
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Figure 3-3: Electrical Circuit Model of Left Ventricle 


The pulmonary circulation and the right ventricle can be modeled similarly. Com- 
bining the four functional units, we get a lumped parameter model shown in Figure 
3-1. The normal values of the resistances and the capacitances in Figure 3-1 can be 
estimated by the relative resistance and volume in a normal person’s cardiovascular 


system. 


39 


3.3. Baroreflex Model 


In order to generate a more realistic simulation, the CVsim adds a baroreflex model to 
simulate the major cardiovascular feedback system. The primary goal of the central 
regulatory mechanisms is to maintain arterial blood pressure P,4 within a very narrow 
range. The arterial baroreflex is the principal mediator of short-term control. The 
control system preserves P, by feeding back transduced values of P4. A feedback 


and control diagram is shown in Figure 3-4. 


vagus Ie 
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Figure 3-4: Feedback and Control Diagram of Baroreflex System. The autonomic 
nervous system (ANS) receives pressure difference then trigger the adjustment of the 
hemodynamic state. 


In our research, since we always assume that we deal with the data when patients 


are in stable state, not in transient state, and thus we may turn off the baroreflex. 


3.4 Some Example Cases 
By using CVsim, we can simulate the hemodynamic state of different ICU patients: 


e The primary effect of mitral stenosis is increased R;;. The compensatory effects 
of mitral stenosis can be simulated by increasing heart rate, increasing total 


blood volume, and increasing the contractility of the left and right ventricles. 
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e The primary effect of vasodilation is decreased R,. The compensatory effects 
can be simulated by increasing heart rate, increasing total blood volume, and 


increasing the contractility of the left ventricle. 


e The primary effect of hypertension is increased R,. The compensatory effects 
can be simulated by increasing heart rate, decreasing venous zero-pressure vol- 
ume, decreasing left ventricular diastolic capacitance, decreasing left ventricular 
systolic capacitance, decreasing venous zero-pressure volume, and increasing ar- 


terial resistance. 


e The primary effect of massive bleeding is decreased total blood volume. The 
compensatory effects can be simulated by increasing heart rate, and decreasing 


left ventricular systolic capacitance. 


e The primary effect of septic shock is decreased R,. The compensatory effects 
can be simulated by increasing heart rate, increasing total blood volume, and 


increasing the contractility of the left ventricle. 


e The primary effect of tricuspid stenosis is increased R,. The compensatory 
effects can be simulated by increasing heart rate, and increasing total blood 


volume. 


In the following chapters, we will develop a method to automatically identify 
the actual parameters of the model based on measurable data such as arterial blood 
pressure, central venous pressure, cardiac output, etc. Due to the limitation of the 
configuration of the circuit, there are some patient cases which are not possible to sim- 
ulate, such as mitral regurgitation, tricuspid regurgitation, ventricular septal defect, 


etc. 
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Chapter 4 


Reasoning Model 


In this chapter, we will talk about the methodology used in our expert system. Tra- 
ditionally, the first issue coming out of building an expert system is how to organize 
and represent knowledge. That is a typical scenario when people used to give qual- 
itative results. In recent years, with the rapid development of mathematical tools, 
many researchers started to use pure or mainly pure numerical methods to try to get 
quantitative results. During the investigation of our current system, we tried unsuc- 
cessfully to use numerical methods, such as the homotopy method, to quantitatively 
solve the problem. Eventually, a new non-numerical method was proposed based on 


the analysis of medical decision making process. 


4.1 A Typical Expert System 


The process of building an expert system is often called knowledge engineering. It 
typically involves a special form of interaction between the knowledge engineer and 
experts in specific domains. The knowledge engineer “translates” the experts’ knowl- 
edge, procedures, and strategies for problem solving into a computer program, as 
shown in Figure 4-1. 

The ideal result will be a computer program successfully solving problems in the 
same manner as human experts. The heart of an expert system is the corpus of knowl- 


edge that accumulates during system building. In expert system design, the main role 
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Figure 4-1: Knowledge Engineering, from [38] 


players are the expert system, the domain expert, the AI consultants, the knowledge 
engineer, the expert system building tools, and the user. Their relationships to each 
other are shown in Figure 4-2. 

The domain expert is a knowledgeable person with a reputation for producing right 
results in a specific domain. The knowledge engineer is a person with a background 
in Al, and is in charge of collaborating with domain expert and AI consultants. The 
AI consultants are respected persons in AI domain, who can always provide guidance 
to the knowledge engineer regarding general expert system questions. In my research, 
Prof. Roger Mark has been serving as the domain expert in cardiovascular diseases, 
Prof. Peter Szolovits has provided some interesting thoughts in searching for the 
methodology for the system as the AI consultant. I am the knowledge engineer. The 


system platform is the ANSI C programming language on a normal UNIX system. 


4.1.1 Representation of Knowledge 


The heart of an expert system is its corpus of knowledge, as we mentioned before. 
The three most widely used knowledge representation in current expert systems are 


rules, semantic nets, and frames. 
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Figure 4-2: Roles in Expert System Building, adapted from [38] 


e A rule-based representation uses the common statement JF... THEN .... 
The matching of rule IF portions to the facts can produce inference actions. 
Rules provide a natural way to describe problem-solving strategies. We can 
think about a scenario when we see a doctor and tell him what the problem is. 
For example, we might state: 

“T have abdominal pain ..., I have ...” 

The doctor will make the diagnosis based on our statements and some additional 
lab tests. How do doctors make their diagnosis? One way to think about is 
that they usually come up with some reasoning processes like: 

“ .. because he has abdominal pain ... he may have stomach problems ...” 
Actually, the physicians’ inference actions result from a rule such as: 


“Tf there is abdominal pain, then there is a possibility for the patient to have 
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stomach problem.” 

This reasoning technique is called forward chaining, because the system uses 
information on the left-hand side (JF part) to generate new information on the 
right side (then part). In the situation of a very large expert system with many 
rules, sometimes, backward chaining is more cost effective. With this inference 
method, the system starts with what it wants to verify and only executes rules 


that are relevant to proving it. 


e A semantic-net-based representation uses a semantic net to represent knowl- 
edge based on a network structure. A semantic net consists of nodes connected 
by arcs that describe the relations between nodes. The isa and has-part relations 
are basic relations in a semantic net. Semantic nets are very useful to represent 


knowledge in domains using taxonomies [38]. 


e A frame-based representation uses a frame to represent knowledge. A frame 
is very much like a semantic net. In a frame, the topmost nodes represent general 
concepts, and the lower nodes represent some instances of the concepts. The 
concept at each node is defined as a collection of attributes and slots that are 
values of those attributes. Each slot can have several procedures attached to it. 
There are three commonly used procedures: I[f-added, If-removed, and If-needed 


procedures [38]. 


4.2 Problem Statement 


In our present study, we used the CV simulator to generate “pseudo-clinical” test 
data which was representative of a variety of disease states. This approach simplified 
the task of designing the expert system, and also provided a quantitative method to 
evaluate the performance of our expert system. Since the actual input parameters of 
the CV model are known for each test case, they can be compared directly with those 
derived from the expert system. 


The CV simulator is initialized by specifying each of the 23 model parameters. 
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The simulator then generates the resultant pressures, flows, and volumes at all sites in 
the CV system. The raw waveforms are then pre-processed by the “feature-detector” 
which derives a set of 13 clinically observable parameters (features) such as heart 
rate, mean arterial blood pressure, pulse pressure, mean central venous pressure, 
pulmonary arterial wedge pressure, cardiac output, etc. (Note that clinically non- 
observable parameters are not included in the feature set.) This feature set is used to 
characterize the physiologic data generated by the model. The diagram of generating 


and pre-processing data is shown in Figure 4-3. 


Pere ee ee te eo eee rer Renee Simulated Waveforms 
-———_1——_ > 
Controlled CV ; Feature 
Parameters Simulator Ls i SCDetector 
>| Data fed to 
the expert system 


Sec cecccstece Sse cececoec cece cece eee eee eee 


Figure 4-3: Diagram of Generating and Pre-processing Data 


The input to the expert system is the 13 dimensional feature set from the patient 
or “pseudo-patient”. Now the task of our expert system has been defined as how 
to derive a set of 23 control parameters for the CV model which result in output 


waveforms and features which match those of the patient. 


4.3 A Pure Numerical Method for Quantitative 
Assessment 


Since our task is to derive 23 control parameters for the CV model, we can think of 
the problem as a 23 dimensional search problem. This led us to think about applying 
some traditional pattern recognition methods to this problem. 

In [40], in order to solve some non-linear algebraic equations, the homotopy 


method was used. When we started to look for a good numerical method to do 
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this 23 dimensional search, we applied the homotopy method to our problem. In the 
following sections, I will start with an introduction to the homotopy method, then 
discuss how we applied this method to our problem, and some lessons learned from 


the work. 


4.3.1 Introduction to the Homotopy Method 


The concept of homotopy was originally introduced in topology. It is a very useful 
concept in categorizing different geometric shapes. Then this method was used in 
computer vision research, typically in object recognition. Morgan et al. in [26] first 
proposed to extend the homotopy concept into solving non-linear algebraic equations. 
In algebra, they define a homotopy as a schedule for transforming the start system (or 
equations) into the target system (equations). The homotopy process can be written 


as follows: 


AH(y,t) = (1—1)G(y) + tF(y) (4.1) 


where G(y) is the starting system and F'(y) is the target system. Here y represents 
a set of variables. Suppose we know the solutions of G(y) (we can always assume so 
because we can make up some algebraic equations with known roots) and we want to 
get the solutions of F'(y). Note that at t = 0, the solutions of H(y,t) = 0 are those 
of G(y) = 0, and at t = 1 the solutions of H(y,t) = 0 are those of F(y) = 0. The 
central idea of homotopy method is if we can track the changes of the solutions from 


H(y,0) to H(y,1), we will be able to solve F(y), as shown in Figure 4-4. 


Original Equations 
H(y,0) = Gly) 


Known Solutions 


Target Equations 
H(y,1) = Fly) 


Unknown Solutions 
9 


HOMOTOPY 


Figure 4-4: Homotopy Method for Solving Non-linear Equations 


AT 


Using Eq. 4.1 to track a path from a known solution (y°,t°), we first predict the 
solution for t = t° + At and then correct the prediction using Newton’s method with 
t fixed. For small Ay and At, the Taylor series of expansion for H gives: 

A(y + Ay,t + At) = A(y,t) + JyAyt+ At (4.2) 

where J, and J; are the Jacobians of H with respect to y and ¢. In the prediction 
step, we have: 

A(y°,t?) =0 (4.3) 


and 


H(y° + Ay, t® + At) =0 (4.4) 


Then we can plug Eq. 4.3 and Eq. 4.4 into Eq. 4.2 and get: 


Ay=—J,*F,At (4.5) 


Since Eq. 4.5 is only an approximate solution, we correct the solution at the new 


value of t by setting At = 0 to get a correction step: 


Ay = —J,'H(y,t) (4.6) 


We can repeat the correction step several times till we feel comfortable with the 
maximum error. Then we can incrementally increase t again, and go back to Eq. 4.2 


and calculate the solutions under new t. 


4.3.2 Howto Apply the Homotopy Method to Our Problem 


As we have discussed before, we can consider our problem as a 23 dimensional pattern 
search problem. By using the homotopy method, we can construct a system with 


known input parameters. For example, we can always start from the normal CV 
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simulator parameters. What we know is the output feature of patients or pseudo- 
patients. We can incrementally make an effort to closely match the output features, 
and meanwhile, track the changes of the CVsim control parameters. Ideally, we will 


get the patient parameters setting at the end. The diagram is shown in Figure 4-5. 
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Figure 4-5: Diagram of Homotopy Method for Searching the Pattern of Patients 


In our system, suppose 7@ = (21,2%2,---,2n) is a vector representing the input 
systemic parameters of CVsim and 7 = (y1, y2,---, yn) is a vector representing the 


output features. The procedures of applying homotopy method is as follows: 


1. Set a new yt, anticipate the new (Az,, Avg,---, Avy )) by solving the linear 


equations: 
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3. Run CVsim and get y! = (445 Yos a a: 
4. Set Ay! = y! — y@. 


5. Refine new 2’ by solving the linear equations: 


Ay, © ee aAry 
0 
Ay, © srlaAer + gelaAt, +--+ airlaAen (4.8) 


6. Repeat step 2 to step 5 until satisfied. 


7. Check whether yf has reached the target yf, if yes, then stop; otherwise go back 
to step 1. 


One of the problems is that we do not have the oe ~ function directly available. 
We need to run CVsim to get a huge 23 dimensional eta table to describe the 


partial derivative function. 


4.3.3. Some Lessons Learned 


We implemented the homotopy method for solving the problem. Unfortunately, it 
did not work as well as we expected. We have learned some lessons from it. There are 
several concerns with using a pure numerical method on the model-based reasoning 


system: 


e The method works very well if only 1 or 2 systemic parameters have been 
changed. This problem seems a traditional AI problem (System only works 
for the simple cases). But the issue here is not that simple. The basic idea 
of the homotopy method is more or less like a gradient method. What is the 
disadvantage of using a gradient method to solve some non-linear problems? As 
we all remember, the major concern is the local minimum vs. global minimum. 


Of course, we can try to add some techniques, such as simulated annealing, to 


a0) 


try to avoid the local minimum. Stimulated annealing allows the search to take 
some uphill steps to escape from the local minimum. Theoretically, system will 
escape from a local minimum, maybe after many steps. In realistic situations, 


those techniques are not applicable because of time issues and stability issues. 


e By using this numerical method, we must ask ourselves a question: how much 
can we trust the CVsim? The reason why this comes out as an issue is that in 
order to do the gradient based analysis, we need to get a huge partial derivative 
look-up table. This look-up table can only be generated from running the 


CVsim. 


e The numerical method is always trying to minimize the total error. However, 


at the same time, it ignores two main problems: 


— For each reasoning (moving) process, it does not follow some common 
sense strategies. In other words, the method does not consider the real 


physiological status. This is why it may give out some ridiculous answers. 


— As we will find out in the next chapter, the decrease of the total error of 
the output features between two systems does not necessarily mean that 


the input parameters are getting closer. 


4.4 Design of the Current System 


All is not lost with the failure of the homotopy method. It reminded us to look for 
some less numerical and more heuristic methods to deal with the reasoning process. 
Our goal is to produce quantitative interpretation or at least semi-quantitative inter- 
pretation. Why is the quantitative analysis much more difficult than the qualitative 
analysis? If I make an analogy to the problem of searching for an apple in a room, 
quantitative analysis means that we need to be able to figure out not just the direc- 
tion of the move, but also how much the move should be. Intuitively, it is almost 
impossible to reach the apple by a single step move (we all have to move and then 


look, then move, ...), unless we are very lucky. So in order to perform quantitative 
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analysis, we can get some feeling about the design strategies from grabbing an apple 
in aroom. That is, the method must be iterative, must be able to tell me something 
at each step, which is related to how far away toward our target after each move. In 
other words, our expert system must be able to run CVsim iteratively, and search for 
the best match. Then the diagram of the current system was put forward in [41], as 


shown in Figure 4-6. 
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Initial Parameters 
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Figure 4-6: Diagram of Our Expert System 


We first establish an initial value for each of the systemic input parameters. Then 
we run CVsim. From the output variables of CVsim, such as mean arterial blood 


pressure, mean central venous pressure, systemic arterial pulse pressure, we use the 
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inference engine to determine which input parameter(s) need to be changed in order 
to decrease the differences between the outputs from CVsim and the pseudo-patient. 
Again, after adjusting the input systemic parameters, we re-run CVsim. This iteration 
process is continued until we feel satisfied with the difference between the outputs 
from CVsim and the pseudo-patient. The comparison between them is made by using 


the following error function: 


__ |Pseudopatient Output Value — CV sim Output Value| 


E, 
‘ Pseudopatient Output Value 


For example, the error function of the mean CVP (Central Venous Pressure) is 
defined as follows: 


|Pseudopatient mean CV P — CV sim mean CV P| 


Ee = ; 
y Pseudopatient mean CV P Value 


In fact, if we look at the Figure 4-6 and compare with the strategy we developed 
in the homotopy method, we may notice that we actually inherited the concept of 
homotopy. The only major difference is that we use heuristic methods to design the 


inference engine instead of pure mathematical methods. 


4.4.1 Initial Values of Systemic Parameters 


In our system, we always start with calculating some computable systemic parameters, 


including: 
1. Heart Rate (beats/min). 


2. Systemic Microvascular Resistance R, (mmHg*sec/ml): 


ie Mean SABP — Mean CVP 
7° Cardiac Output /60 


ABP: Systemic Arterial Blood Pressure; CVP: Central Venous Pressure. 
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3. Pulmonary Microvascular Resistance Rya (mmHg*sec/ml): 


Pics Mean PAP — Mean PCWP 
= Cardiac Output /60 


PAP: Pulmonary Arterial Pressure; PCWP: Pulmonary Capillary Wedge Pres- 


sure. 
4. Systemic Arterial Capacitance C, (ml/mmHg): 


Cs Stroke Volume 
“~~ Pulse Pressure 


Note that we do not calculate some resistances that look calculable, such as Right 


Ventricular Outflow Resistance. Theoretically, we can use 


_ Mean RV P — Mean PAP (when RVP > PAP) 


ie Cardiac Output /60 


to estimate R,, (PAP: Pulmonary Arterial Pressure, RVP: Right Ventricular Pres- 
sure). The problem with the above formula is that the normal value of this resistance 
is very small, equal to 0.003PRU. Assume that this person’s cardiac output is nor- 
mal (about 5000cc/min), then the mean pressure difference during systole between 


pulmonary arterial pressure and right ventricular pressure will be: 


5000(ml/min) 


: el 1) * —————_—— 
0.003(mmHg * sec/ml) x 0(sec/min) 


= 0.25mmAg 


We can not assume that in real clinical areas, 0.25mmHg is detectable. So unless 
there is severe valvular stenosis symptom, I assume that these resistances are normal. 
As a matter of fact, physicians usually make the decision about pulmonary valve 
stenosis based on the instantaneous pressure gradient between the right ventricle and 
the pulmonary artery. 

There are also several parameters which do not change based on our assumptions. 


They are: 
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1. Transthoracic Pressure: The changes of the transthoracic pressure can change 
the relative filling pressure. Ignoring the changes of transthoracic pressure, we 
eliminate the effects of the respiratory system. So our system can not be used 


for the patients on ventilators. 


2. All Zero-pressure Volumes: The changes of the zero-pressure volumes can be 
compensated by changing the total blood volume, because from the hemody- 
namic point of view, we are concerned with the effective blood volume (=total 


blood volume - the summation of zero-pressure volumes). 


3. Venous Capacitance: Clinically, incremental venous capacitance does not change 


in most of the patient cases. 


All of the remaining input CVsim parameters are initialized to normal default 


values. 


4.4.2 Flow Chart of the Reasoning Process 


In our system, medical knowledge is represented as rules. We have two flow charts in 


our system: 


e Figure 4-7, Figure 4-8 and Figure 4-9 show the reasoning process for patients 


with no valvular disease. 


e Figure 4-10, Figure 4-11 and Figure 4-12 show the reasoning process for ana- 


lyzing patients with possible mitral stenosis and/or tricuspid stenosis diseases. 


To determine which flow chart the system uses for each patient, we start to run 


our expert system while asking some questions, such as: 


e Is there any diastolic murmur? 


e Is there any sign of mitral stenosis from the echocardiogram? 
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The detection of tricuspid stenosis is much easier than the detection of the mitral 
stenosis, since we have the measurement of right ventricular pressure. We can make 
the decision based on the instantaneous pressure gradient between the central venous 


pressure and the right ventricular pressure. 


4.4.2.1 Flow Chart for Non-valvular Disease Patients 


Figure 4-7 shows the top level of the reasoning process. It always starts from checking 
the difference of the cardiac output between the pseudo-patient and the CVsim. The 


comparison gives three possible results: 


e High means the cardiac output of the CVsim is higher than that of the pseudo- 


patient. 


e Equivalent means the cardiac output of the CVsim is within the same range 


(5% for cardiac output) with that of the pseudo-patient. 


e Low means the cardiac output of the CVsim is lower than that of the pseudo- 


patient. 


Since we assume there is no valvular disease in this flow chart, there is no action to 
change the input and output flow resistances. The possible changing systemic param- 
eters are: total blood volume, left ventricular systolic capacitance, right ventricular 
systolic capacitance, left ventricular diastolic capacitance, right ventricular diastolic 
capacitance, systemic arterial capacitance and pulmonary arterial capacitance. In 
Figure 4-7, tune up means the situation of the pseudo-patient’s heart is better than 
that of the CVsim. In this case, we need to tune up the systemic parameters of CVsim. 
Tune down is the opposite situation of tune up. Figure 4-8 continues the reasoning 
process of tune up or down. In Figure 4-7, adjust means we need to change both the 
systolic capacitance and the diastolic capacitance of left/right ventricle. Figure 4-9 
continues the reasoning process of adjustment. 

Finally, if all the measurable parameters of CVsim are considered to be equivalent 
with that of the pseudo-patient, the iteration should stop. It means the current 


CVsim’s systemic parameters can be used to simulate the status of the pseudo-patient. 
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4.4.2.2 Flow Chart for Patients with Possible Mitral Stenosis and Tricus- 


pid Stenosis Diseases 


From the echocardiogram, we may find some evidence of mitral stenosis and/or tricus- 
pid stenosis. The possible changing systemic parameters are: total blood volume, left 
ventricular systolic capacitance, right ventricular systolic capacitance, left ventricular 
diastolic capacitance or inflow resistance, right ventricular diastolic capacitance or 
inflow resistance, systemic arterial capacitance and pulmonary arterial capacitance. 
Note that in the case of mitral stenosis, pulmonary wedge pressure is not equal to the 
left ventricular end diastolic pressure. Figure 4-11 continues the reasoning process of 
tune up or down. Figure 4-12 continues the reasoning process of adjustment. If there 
are signs of mitral stenosis, we change the value of left ventricular inflow resistance 
instead of the value of the left ventricular diastolic capacitance. If there are signs of 
tricuspid stenosis, we change the value of right ventricular inflow resistance instead 


of the value of the right ventricular diastolic capacitance. 


4.4.3 Some Heuristic Components in Our Expert System 


In designing an expert system, there are a number of other important factors that 
we need to think about besides which parameter(s) need to be changed. A very 
important factor in designing the search control process is how to define the step size 
for each legitimate move. An easy approach is to set a very small fixed step size. The 
problem with the small fixed step size is obvious: it may take us a very long time to 
get to the target point from the original point. So in our system, we decided to use 
variable step sizes instead. The way that we set/change step size during the search 


control process is described as follows: 


e We set initial unit step sizes, which usually tend to be conservative, when we 


start the search control process. 


e The step size is proportional to the difference between the outputs of the CVsim 


and the pseudo-patient. For example, we change the systemic arterial capaci- 
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tance based on the difference of the pulse pressure of the CVsim and the pseudo- 


patient. The larger the difference is, the larger is the step size. 


Step Size of Systemic Arterial Capacitance = 
|Pseudopatient Pulse Pressure — CV sim Pulse Pressure| * 


Unit Step Size 


e We set an increase unit step size and a decrease unit step size separately for each 
parameter. The reason for this strategy is that we want to avoid oscillation. The 
disadvantage of defining only one unit step size for both decrease and increase 


is shown in Figure 4-13. 


e If the new proposed change for a parameter is exactly in the same direction as 
the last change, we increase the appropriate unit step size by 50%. This is what 


we called “Reward” strategy, as shown in Figure 4-14. 


e If the new proposed change for a parameter is in the opposite direction with 
the last change, we will go back to the place where we were, then decrease the 
appropriate unit step size by 50%. This is what we called “Punish” strategy. 


Figure 4-15 shows how it works. 
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Figure 4-7: Top Level Flow Chart of Reasoning Process for Non-valvular Dis- 
ease Patients. C.O.=Cardiac Output; LVEDP=Left Ventricular End-Diastolic Pres- 
sure; ABP=Mean Systemic Arterial Blood Pressure; LEF=Left Ejection Frac- 
tion; CVP=Mean Central Venous Pressure; PAP=Mean Pulmonary Arterial Pres- 
sure; REF=Right Ejection Fraction; SAPP=Systemic Arterial Pulse Pressure; 
PAPP=Pulmonary Arterial Pulse Pressure; Equv.=Equivalent. Continued with Fig- 
ure 4-8 and Figure 4-9. 
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or down? 
High Low High Low 
Decrease Increase Increase Decrease 
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Figure 4-8: Tune Up/down Reasoning Process for Non-valvular Disease Patients. 
Continued for Figure 4-7 when we need to tune or down the ventricular functions. 
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Figure 4-9: Adjust Reasoning Process for Non-valvular Disease Patients. Continued 
for Figure 4-7 when we need to adjust the ventricular capacitances based on ejection 
fraction. 
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Figure 4-10: Top Level Flow Chart of Reasoning Process for Patients with Possi- 
ble Mitral Stenosis and/or Tricuspid Stenosis. C.O.=Cardiac Output; LVEDP=Left 
Ventricular End-Diastolic Pressure; PCWP=Pulmonary Capillary Wedge Pres- 
sure; ABP=Mean Systemic Arterial Blood Pressure; LEF=Left Ejection Frac- 
tion; CVP=Mean Central Venous Pressure; PAP=Mean Pulmonary Arterial Pres- 
sure; REF=Right Ejection Fraction; SAPP=Systemic Arterial Pulse Pressure; 
PAPP=Pulmonary Arterial Pulse Pressure. Continued with Figure 4-11 and Figure 
4-12. 
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Figure 4-11: Tune Up/down Reasoning Process for Patients with Possible Mitral 
Stenosis and/or Tricuspid Stenosis. Continued for Figure 4-10 when we need to tune 
or down the ventricular functions. 
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Figure 4-12: Adjust Reasoning Process for Patients with Possible Mitral Stenosis 
and/or Tricuspid Stenosis. Continued for Figure 4-10 when we need to adjust the 
ventricular capacitances based on ejection fraction. 
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Parameter Value 


Correct Value 


step size 


Number of Iterations 
Figure 4-13: Why do we need two unit step sizes for each parameter? As we can see, 


we may overshoot sometimes. If there is only one unit step size for both increase and 
decrease, that may cause the system to oscillate. 


Parameter Value 


1 2 3 4 Number of Iterations 


Figure 4-14: Reward Strategy. If the proposed change is in the same direction as the 
last change, we tend to think that we are still quite far away from our target, we need 
to be a little bit more aggressive in order to save iteration times. 
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Figure 4-15: Punish Strategy. It is another strategy to avoid oscillation. 
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Chapter 5 


Development and Evaluation of 


the System 


Evaluating an expert system is very difficult because sometimes, there may be no 
formal way to prove a given answer is correct or acceptable or the best possible. The 
evaluation process usually involves consulting other experts (not the one where we 
had most of the domain knowledge) in the domain. Ultimately, user acceptance can 
make or kill an expert system. Obviously, users will not accept an expert system with 
poor performance; neither will they accept a high performance expert system that 
gives confusing answers, or is hard to learn how to use, or takes a long time to get 
results. 

A number of evaluations of medical expert systems have been put forward in 
[19][22][25], most of them were focusing on evaluation of the qualitative analysis. In 
these evaluation processes, the results from expert systems are presented to human 
experts, and they are asked to fill a sheet for the evaluation for each case. In [22], 


reviewers were given the following 5 choices: 


e Correct means the result is a reasonable accounting of the findings. 


e Possible means the result is a possible interpretation of the findings, but not 


the best. 
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e Partly Correct means the result is mostly reasonable, but has some minor 


mistakes. 
e Wrong means the result is wrong. 


e Seriously Wrong means the result is wrong and dangerous. 


It is even more difficult to evaluate the performance of our expert system, because 
our system provides quantitative results. An ideal evaluation of our program would 
begin by identifying a large number of actual patient cases that occur in the cardio- 
vascular domain covered by this program. In this approach to evaluation, the actual 
patient data play two roles. First, they tell us how frequently various conditions arise, 
and therefore how important the program’s performance is on such cases. Second, 
they provide us a natural variety of distinct cases against which to test our program. 

Performing such an evaluation is quite difficult because of the need for a large 
number of cases. In addition, we believe that relatively small variations in the pseudo- 
patient case would not yield significantly better or worse performance by our program. 
Therefore, instead of testing on cases derived from actual patient data, we have eval- 
uated our approach by simulating a varied set of clinically important and distinguish- 
able cases that span the domain of expertise of the program. These cases can then 
illustrate the strengths and weaknesses of the program, although we cannot expect 
them to yield a statistical estimate of how frequently it would do well or poorly with 
situations abstracted from actual cases. Of course, such illustrations are exactly what 
we need in a formative evaluation, whose major purpose is to identify areas in which 
additional work will be needed. 

While we were developing the reasoning and knowledge base of our system, we 
used a learning set of 16 simulated cases for both testing and refining the algorithms. 
In the next section, we will present and discuss the test results of our expert system 


on these cases. 
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5.1 Results and Discussion 


All of the relevant tables and figures can be found in Appendix A. In this section, we 
will analyze the results and discuss potential improvements for our system. For each 


case, we will present two tables: 


e Comparison of Pseudo-Patient Input Parameters with Estimated Input Param- 


eters 
e Comparison of Output Variables of Pseudo-Patient with those of Simulator 
Two figures will also be provided for each case: 


e Total Output Error Function vs. Number of Iterations: The total output error 


function is defined as: 


Eourpur = VE? +-::+ Bb? +-+++ EF? 


_ |Pseudopatient Output Value — CV sim Output Value| 


E, 
: Pseudopatient Value 


e Total Input Error Function vs. Number of Iterations: The total input error 


function is defined as: 


Ernwpur = VE? +-++++ E24+-+++ FE? 


E= |Pseudopatient Parameter Value — Estimated Parameter Value| 
Pseudopatient Parameter Value 


Note that we have 7 volume related input parameters for CVsim, but since we can 
only characterize the system state based on the total effective blood volume (=total 
blood volume - summation of all zero-pressure volumes), only the input error of the 
total effective blood volume is taken into account for calculating the total input error. 

In addition, it is also very important to conduct a qualitative evaluation of the 


performance of the system. The domain expert (Dr. Roger Mark) who had helped 
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to design the system was asked to participate in the qualitative evaluation. The 
system’s estimated input parameters and CVsim’s output variables were presented 
to Dr. Mark, and he was asked to provide a corresponding clinical diagnosis. The 


expert system’s performance was then rated using the following scale: 


e Excellent means that: 
1. The diagnosis implied by the expert system’s estimated parameters is con- 
sistent with the original diagnosis modeled by the test case. 
2. The algorithm converges. 


3. The final total input error is less than 20%. 
e Acceptable means that: 


1. The diagnosis implied by the expert system’s estimated parameters is con- 
sistent with the original diagnosis modeled by the test case. 
2. The algorithm does not converge and/or the (final)! total input error is 


greater than 20%. 


e Unsatisfactory means that the diagnosis implied by the expert system’s esti- 
mated parameters is not consistent with the original diagnosis modeled by the 


test case, and would probably lead to an erroneous diagnosis in clinical areas. 


5.1.1 Cases without Compensatory Effects 


We start our discussion with the testing results on the cases without compensatory 
effects. These cases represent relatively simple physiological changes with respect to 


normal people. 


If the algorithm does not converge, then there is no final total input error. 
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5.1.1.1 Blood Loss 


This case represents a patient with blood loss. The status of a blood loss patient was 
simulated on CVsim by decreasing total blood volume. After running about 4 minutes 
(8 iterations) on a Sun SPARC-10 station, the iteration stops. Table A.1 shows the 
comparison of the input parameters of the pseudo-patient to those estimated by the 
expert system, with reference to their normal values. Table A.2 shows the comparison 
of the output variables of the pseudo-patient to those of the simulator, with reference 
of their normal values. Figure A-2 shows the output error function for each iteration. 
Since we know how we set the input parameters of CVsim, we also printed out how the 
input error function behaves for each iteration in Figure A-1. The final total input 
error is 7.4%, and the final total output error is 12.0%. In this case, the algorithm 
converges very well. The rating of this case given by the reviewer is excellent. 

By comparing Figure A-2 and Figure A-1, we may notice that the drop of the 
output error function does not necessarily correspond to the drop of the input error 
function. For example, from the first iteration to the second iteration, there is a huge 
drop in the output error function; however, the input error function increases. This 
may be part of the reason why the homotopy method failed to solve this problem, 
because the basis of the homotopy method is how to achieve a smaller output error 
function for each step. In our system, we use total error function to determine whether 
we should stop the iteration, but we do not use the output error function to determine 


whether the move for each iteration is legitimate or not. 


5.1.1.2 Arteriolar Dilation 


This case represents a patient with arteriolar dilation, and was simulated on CVsim 
by decreasing peripheral resistance. After running about 2 minutes (4 iterations) on 
a Sun SPARC-10 station, the iteration stops. Table A.3 shows the comparison of 
the input parameters of the pseudo-patient to those estimated by the expert system, 
with reference to their normal values. Table A.4 shows the comparison of the output 


variables of the pseudo-patient to those of the simulator, with reference to their 
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normal values. Figure A-4 shows the output error function for each iteration. The 
input error function for each iteration is shown in Figure A-3. The final total input 
error is 1.3%, and the final total output error is 1.6%. In this case, the algorithm 


converges very well. The rating of this case given by the reviewer is excellent. 


5.1.1.3 Left Ventricular Systolic Dysfunction 


This case represents a patient with left ventricular systolic dysfunction. The status 
of this patient was simulated on CVsim by increasing the left ventricular systolic 
capacitance. After running about 3 minutes (7 iterations) on a Sun SPARC-10 station, 
the iteration stops. Table A.5 shows the comparison of the input parameters of the 
pseudo-patient to those estimated by the expert system, with reference to their normal 
values. Table A.6 shows the comparison of the output variables of the pseudo-patient 
to those of the simulator, with reference to their normal values. Figure A-6 shows the 
output error function for each iteration. The input error function for each iteration is 
shown in Figure A-5. The final total input error is 3.5%, and the final total output 
error is 2.8%. In this case, the algorithm converges very well. The rating of this case 


given by the reviewer is excellent. 


5.1.1.4 Arteriolar and Venous Dilation 


This case represents a patient with arteriolar and venous dilation. The case was 
simulated on CVsim by increasing venous zero-pressure volume, and decreasing pe- 
ripheral resistance. After running about 4 minutes (9 iterations) on a Sun SPARC-10 
station, the iteration stops. Table A.7 shows the comparison of the input parame- 
ters of the pseudo-patient to those estimated by the expert system, with reference to 
their normal values. Table A.8 shows the comparison of the output variables of the 
pseudo-patient to those of the simulator, with reference to their normal values. Figure 
A-8 shows the output error function for each iteration. The input error function for 
each iteration is shown in Figure A-7. The final total input error is 13%, and the 
final total output error is 5.5%. In this case, the algorithm converges very well. The 


rating of this case given by the reviewer is excellent. 
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5.1.2 Cases with Compensatory Effects 


The following cases represent some pathological cases with compensatory effects. In 
these cases, primary changes are caused by the disease directly, and the compensatory 


effects are caused by various neural and humoral compensatory mechanisms. 


5.1.2.1 Hypertension 


This case represents a patient with hypertension. The primary effect was simulated 
on CVsim by increasing peripheral resistance. The compensatory effects were in- 
creased heart rate, decreased venous zero-pressure volume, increased left ventricular 
contractility, decreased left ventricular diastolic compliance, and decreased arterial 
compliance. After running about 10 minutes (24 iterations) on a Sun SPARC-10 
station, the iteration stops. Table A.9 shows the comparison of the input parame- 
ters of the pseudo-patient to those estimated by the expert system, with reference to 
their normal values. Table A.10 shows the comparison of the output variables of the 
pseudo-patient to those of the simulator, with reference to their normal values. Figure 
A-10 shows the output error function for each iteration. The input error function for 
each iteration is shown in Figure A-9. The final total input error is 10.0%, and the 
final total output error is 9.7%. In this case, the algorithm converges very well. The 


rating of this case given by the reviewer is excellent. 


5.1.2.2. Right Ventricular Infarction 


This case represents a patient with right ventricular infarction. The primary effect 
was simulated by increasing right ventricular systolic capacitance and decreasing right 
ventricular diastolic capacitance. The compensatory effects were increased heart rate 
and increased peripheral resistance. The program did not stop by itself because it 
stayed in a local minimum state. Table A.11 shows the comparison of the input 
parameters of the pseudo-patient to those estimated by the expert system during the 
24th iteration (system takes about 10 minutes to run 24 iterations). Table A.12 


shows the comparison of the output variables of the pseudo-patient to those of the 
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simulator during the 24th iteration, with reference to their normal values. Figure 
A-12 shows the output error function for each iteration. The input error function for 
each iteration is shown in Figure A-11. The algorithm did not converge in this case. 
The rating of this case given by the reviewer is acceptable. 

From Table A.12, we can see that just by comparing the output variables during 
the 24th iteration, the simulator output is very much like the pseudo-patient’s output. 
It means that if we loosen our comparison boundary (with the sacrifice of the precision 
of results), the program could have stopped by itself. From Table A.11, we notice 
that the problem is that our system tended to “think” that the right ventricular 
systolic function is worse than it should be, and the right ventricular diastolic function 
is better than it should be. The system gets stuck in this local minimum state. 
In order to avoid the local minimum state, we may apply a similar strategy with 
simulated annealing [30], which allows us to impose a probabilistic value for each 
branch in the decision tree to increase the uncertainty during each iteration. For 
example, if the mean pulmonary arterial pressure of the CVsim is higher than that 
of the pseudo-patient, the current algorithm will always increase the right ventricular 
systolic capacitance. By imposing a probabilistic value (for example: 80%) to this 
rule, the algorithm may execute the rule by 80% chance or continue the search process 


to check the right ventricular ejection fraction by 20% chance. 


5.1.2.3. Hypertrophic Cardiomyopathy 


This case represents a patient with hypertrophic cardiomyopathy without outlet ob- 
struction. The primary effect was simulated by decreasing left ventricular systolic 
capacitance, decreasing left ventricular diastolic capacitance, and decreasing right 
ventricular diastolic capacitance. The compensatory effects were increased heart 
rate, increased total blood volume, decreased venous zero-pressure volume, and de- 
creased peripheral resistance. After running about 30 minutes (74 iterations) on a 
Sun SPARC-10 station, the iteration stops. Table A.13 shows the comparison of 
the input parameters of the pseudo-patient to those estimated by the expert system, 


with reference to their normal values. Table A.14 shows the comparison of the out- 
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put variables of the pseudo-patient to those of the simulator, with reference to their 
normal values. Figure A-14 shows the output error function for each iteration. The 
input error function for each iteration is shown in Figure A-13. The final total input 
error is 14.6%, and the final total output error is 13.6%. In this case, the algorithm 
converges well, but not very fast. The rating of this case given by the reviewer is 
excellent. 

From Figure A-14, we can see some interesting behavior of the output error 
function during the 5th iteration to 11th iteration. There is a huge rise of the input 
error function during the 7th iteration, and a huge drop during the 8th iteration, then 
another huge rise during the 9th iteration. That is when our punish strategy started 
to work. During the 7th iteration, the reasoning process increased the total blood 
volume by 2983.42ml. Then during the 8th iteration, the reasoning process tended to 
decrease the total blood volume. However, before it really decreased the total blood 
volume, it looked backward and noticed that the problem was last time, the reasoning 
process increased the total blood volume way too much. Then the program went 
back to the position of the 6th iteration, and increased the total blood volume by 
2237.18ml. Again, we see that the changes of the total output error function does not 


necessarily correspond to the changes of the total input error function. 


5.1.2.4 Pulmonary Hypertension 


This case represents a patient who has pulmonary hypertension. The primary effect 
was simulated by increasing pulmonary resistance. The compensatory effects were 
decreased right ventricular systolic capacitance, increased heart rate, increased total 
blood volume, and decreased pulmonary capacitance. The program did not stop by 
itself, because it stayed in a local minimum state. Table A.15 shows the comparison 
of the input parameters of the pseudo-patient to those estimated by the expert system 
during the 24th iteration. Table A.16 shows the comparison of the output variables 
of the pseudo-patient to those of simulator during the 24th iteration, with reference to 
their normal values. Figure A-16 shows the output error function for each iteration. 


The input error function for each iteration is shown in Figure A-15. The rating of 
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this case given by the reviewer is acceptable. 

The major problems of the result are too much effective blood volume and no 
change of pulmonary arterial capacitance. The problem of the algorithm is because 
of the negative effect of our punish strategy. During the 66th iteration, the algorithm 
decreased the left ventricular systolic capacitance. However, on the 67th iteration, 
the algorithm tended to increase the the left ventricular systolic capacitance. Then 
it noticed that the problem may be because last time the left ventricular systolic 
capacitance had been decreased too much. The algorithm decreased the unit step 
size. Then it was still overshoot. Then according to our punish strategy, the unit 
step size needed to be decreased again. Finally, it stayed in that position which 
was always trying to decrease this capacitance, but never succeeded. To solve this 
problem, we may impose a probabilistic value for each branch in the decision tree to 


increase the uncertainty during each iteration. 


5.1.2.5 Large AV Fistula 


This case represents a patient with large AV fistula. The primary effect was simulated 
by decreasing peripheral resistance. The compensatory effects were increased heart 
rate, increased total blood volume, decreased left ventricular systolic capacitance, de- 
creased right ventricular systolic capacitance, decreased venous zero-pressure volume, 
and decreased pulmonary resistance. After running about 10 minutes (23 iterations) 
on a Sun SPARC-10 station, the iteration stops. Table A.17 shows the comparison of 
the input parameters of the pseudo-patient to those estimated by the expert system, 
with reference to their normal values. Table A.18 shows the comparison of the out- 
put variables of the pseudo-patient to those of the simulator, with reference to their 
normal values. Figure A-18 shows the output error function for each iteration. The 
input error function for each iteration is shown in Figure A-17. The final total input 
error is 55.7%, and the final total output error is 6.1%. In this case, the algorithm 
converges well. The rating of this case given by the reviewer is acceptable. 

The major difference between the results and the set-up parameters is that our 


expert system tended to add extra blood volume and worsen the diastolic functions of 
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both ventricles a bit more. From Table A.18, we can see that the output observable 


parameters are very close. 


5.1.2.6 Dilated Cardiomyopathy and Bi-ventricular Congestive Heart Fail- 


ure 


This case represents a patient with dilated cardiomyopathy and bi-ventricular con- 
gestive heart failure. The primary effect was simulated by increasing left ventricu- 
lar systolic capacitance, increasing right ventricular systolic capacitance, decreasing 
left ventricular diastolic capacitance, increasing left ventricular zero-pressure volume, 
and increasing right ventricular zero-pressure volume. The compensatory effects were 
increased heart rate, increased total blood volume, decreased arterial capacitance, 
increased right ventricular diastolic capacitance, decreased pulmonary arterial capac- 
itance, decreased venous zero-pressure volume, and increased pulmonary resistance. 
After running about 35 minutes (82 iterations) on a Sun SPARC-10 station, the it- 
eration stops. Table A.19 shows the comparison of the input parameters of the 
pseudo-patient to those estimated by the expert system, with the reference to their 
normal values. Table A.20 shows the comparison of the output variables of the 
pseudo-patient to those of the simulator, with reference to their normal values. Fig- 
ure A-20 shows the output error function for each iteration. The input error function 
for each iteration is shown in Figure A-19. The final total input error is 52.0%, and 
the final total output error is 16.8%. In this case, the algorithm converges well, but 
not fast enough. From Figure A-20, we can find the similar behavior of the output 
error function as discussed in hypertrophic cardiomyopathy case. The rating of this 


case given by the reviewer is acceptable. 


5.1.2.7 Acute Myocardial Infarction with Bi-ventricular Heart Failure 


This case represents a patient with an acute myocardial infarction and bi-ventricular 
heart failure. The primary effect was simulated by increasing left ventricular sys- 
tolic capacitance, increasing right ventricular systolic capacitance, and decreasing 


left ventricular diastolic capacitance. The compensatory effects were increased heart 
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rate, decreased arterial capacitance, increased right ventricular diastolic capacitance, 
decreased pulmonary arterial capacitance, decreased venous zero-pressure volume, in- 
creased peripheral resistance, and increased right ventricular zero-pressure volume. 
After running about 33 minutes (77 iterations) on a Sun SPARC-10 station, the iter- 
ation stops. Table A.21 shows the comparison of the input parameters of the pseudo- 
patient to those estimated by the expert system, with reference to their normal values. 
Table A.22 shows the comparison of the output variables of the pseudo-patient to 
those of the simulator, with reference to their normal values. Figure A-22 shows the 
output error function for each iteration. The input error function for each iteration is 
shown in Figure A-21. The final total input error is 14.6%, and the final total output 
error is 10.1%. The rating of this case given by the reviewer is excellent. Again, from 
Figure A-22 and Figure A-21, we may notice that sometimes, even if there is a huge 


rise of the output error, there may not be a rise of the input error. 


5.1.2.8 Septic Shock 


This case represents a patient in septic shock. The primary effect was simulated by 
decreasing peripheral resistance, increasing venous zero-pressure volume, and increas- 
ing left ventricular systolic capacitance. The compensatory effect was increased heart 
rate. After running about 22 minutes (51 iterations) on a Sun SPARC-10 station, 
the iteration stops. Table A.23 shows the comparison of the input parameters of 
the pseudo-patient to those estimated by the expert system, with the reference to 
their normal values. Table A.24 shows the comparison of the output variables of 
the pseudo-patient to those of the simulator, with reference to their normal values. 
Figure A-24 shows the output error function for each iteration. The input error 
function for each iteration is shown in Figure A-23. The final total input error is 
14.2%, and the final total output error is 10.7%. In this case, the algorithm converges 


well. The rating of this case given by the reviewer is excellent. 
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5.1.2.9 Pulmonary Hypertension and Right Heart Failure 


This case represents a patient with pulmonary hypertension and right heart failure 
(cor pulmonale). The primary effect was simulated by increasing pulmonary resis- 
tance, and increasing right ventricular systolic capacitance. The compensatory effects 
were increased heart rate, increased total blood volume, decreased peripheral arterial 
capacitance, decreased venous zero-pressure volume, increased peripheral resistance. 
After running about 23 minutes (50 iterations) on a Sun SPARC-10 station, the it- 
eration stops. Table A.25 shows the comparison of the input parameters of the 
pseudo-patient to those estimated by the expert system, with the reference to their 
normal values. Table A.26 shows the comparison of the output variables of the 
pseudo-patient to those of the simulator, with reference to their normal values. Fig- 
ure A-26 shows the output error function for each iteration. The input error function 
for each iteration is shown in Figure A-25. The final total input error is 6.6%, and 
the final total output error is 12.8%. In this case, the algorithm converges very well. 
From Figure A-26, we can see again the punish strategy has worked. The rating of 


this case given by the reviewer is excellent. 


5.1.2.10 Mitral Stenosis I 


This case represents a patient with mitral stenosis. The primary effect was simulated 
by increasing left ventricular inflow resistance. The compensatory effects were in- 
creased heart rate, increased total blood volume, decreased right ventricular systolic 
capacitance and increased pulmonary resistance. In this case study, we intentionally 
indicated to the algorithm that there is no evidence of mitral stenosis and tried to 
see what results we would get. After running about 13 minutes (30 iterations) on 
a Sun SPARC-10 station, the iteration stops. Table A.27 shows the comparison of 
the input parameters of the pseudo-patient to those estimated by the expert system, 
with reference to their normal values. Table A.28 shows the comparison of the out- 
put variables of the pseudo-patient to those of the simulator, with reference to their 


normal values. Figure A-28 shows the output error function for each iteration. The 
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input error function for each iteration is shown in Figure A-27. The final total input 
error is 120.3%, and the final total output error is 7.3%. The rating of this case given 
by the reviewer is unsatisfactory. 

The major difference of our results with set-up parameters are the left ventricular 
inflow resistance and the left ventricular diastolic capacitance. However, clinically, it 
is also very hard to determine the mitral stenosis case without knowing this person’s 
left ventricular pressure. Our expert system explained the problem as a left ventricular 
filling problem, which is correct because with a very high left ventricular diastolic 


pressure, we usually turn our concern to this person’s diastolic problems. 


5.1.2.11 Mitral Stenosis IT 


This case represents a patient with mitral stenosis. The primary effect was simulated 
by increasing left ventricular inflow resistance. The compensatory effects were in- 
creased heart rate, increased total blood volume, decreased right ventricular systolic 
capacitance and increased pulmonary resistance. In this case study, we indicated 
to the algorithm that there is evidence of mitral stenosis. After running about 11 
minutes (27 iterations) on a Sun SPARC-10 station, the iteration stops. Table A.29 
shows the comparison of the input parameters of the pseudo-patient to those esti- 
mated by the expert system, with reference to their normal values. Table <A.30 
shows the comparison of the output variables of the pseudo-patient to those of the 
simulator, with reference to their normal values. Figure A-30 shows the output error 
function for each iteration. The input error function for each iteration is shown in 
Figure A-29. The final total input error is 16.5%, and the final total output error is 
11.1%. The algorithm converges well in this case. The rating of this case given by 
the reviewer is excellent. 

The current algorithm for analyzing the valvular disease is based on the assump- 
tion that if there is any inflow valvular problem, we assume that the diastolic capac- 


itance is normal. 
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5.1.2.12 Tricuspid Stenosis 


This case represents a patient with tricuspid stenosis. The primary effect was sim- 
ulated by increasing right ventricular inflow resistance. The compensatory effects 
were increased heart rate, and increased total blood volume. After running about 
10 minutes (23 iterations) on a Sun SPARC-10 station, the iteration stops. Table 
A.31 shows the comparison of the input parameters of the pseudo-patient to those 
estimated by the expert system, with reference to their normal values. Table A.32 
shows the comparison of the output variables of the pseudo-patient to those of the 
simulator, with reference to their normal values. Figure A-32 shows the output error 
function for each iteration. The input error function for each iteration is shown in 
Figure A-31. The final total input error is 4.3%, and the final total output error is 
4.6%. The algorithm converges very well in this case. The rating of this case given 


by the reviewer is excellent. 


5.2 Summary of the Testing 


We tested our expert system on 16 typical cases, in which 11 cases were rated as 
excellent, 4 cases were rated as acceptable and 1 case was rated as unsatisfactory. 
For the unsatisfactory case, the reason for the wrong diagnosis is lack of knowledge. 
This evaluation provides us with a reasonable picture of the diagnostic performance 
of our expert system. The number of cases in the evaluation was too small to measure 
the performance relative to physicians. In general, the more complicated the case is, 
will the longer time the system take to figure out the problem. There is no serious 
error made by the expert system, but there are some problems, such as time issues, 
and lack of convergence. They imply that further refinement of the reasoning process 
is needed rather than the fundamental changes to the reasoning mechanisms and 
methodology. Currently, we are pleased with the results and believe that the system 


can be improved. 
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Chapter 6 


Conclusion 


In this thesis, we have developed an expert system for quantitatively or semi-quantita- 
tively analyzing ICU hemodynamic data. We have designed a new iterative approach, 
which uses a cardiovascular simulator as a model. The application of this method- 
ology can also be applied outside the area of medicine to solve system identification 
problems. In our system, clinical rules and other available test results (such as an 
echocardiogram) have been used to guide the direction of the iteration in our system. 
The results of our case studies suggested that this approach is promising for solving 


other clinical problems. There are also some unsolved problems in our system: 


e How can we analyze aortic stenosis and pulmonary stenosis? For pulmonary 
stenosis, we can make the quantitative analysis based on the pressure gradient 
between right ventricular pressure and pulmonary arterial pressure. For aortic 
stenosis, it will be very difficult, because left ventricular pressure monitoring is 
not a routine procedure in ICU areas. Without left ventricular pressure data, 
aortic stenosis can possibly be detected by echocardiogram. However, we need 


to be careful to use the quantitative description provided by the echocardiogram. 


e How to further increase the efficiency of the system? Or how to decrease the 
iteration times? Setting more aggressive unit step sizes may be one way, but 


that may also increase the instability of the system. 
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Due to the limitation of the cardiovascular simulator that we used in our research, 


our system can not be used to analyze some cases such as mitral regurgitation. 


6.1 Practical Consideration 


Of course, eventually, we hope to use our system in ICU areas by feeding the system 
with real patient data. Then there will be several major problems that need to be 


solved: 


e Time issue: what is an acceptable time period for the system to converge on an 


answer? 


e As we have shown in two of our cases, sometimes the system may just stay in 
a local minimum state and may iterate forever. Obviously, we cannot count 
on the end-users to stop the program in that case. Then how to avoid this 
situation? We may set a maximum iteration time, once the total iteration time 
reaches the maximum iteration time, we can just stop the program and choose 
the input parameters which correspond to the minimum output error as our 


estimation. 


e We need to pre-process the real ICU data to get rid of noise, artifacts, etc. 


6.2 Future Related Work 


There are a number of improvements that can be made to our current system: 


e For realizing the real time analysis, we need to have some temporal rules instead 
of static rules. Actually, real time analysis may be easier than figuring out the 
problem for the first time, because we do not expect a dramatic change from 
one minute to another minute. It means that the distance from the “target” to 


the “original” state is closer. 
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e During the period of writing this thesis, I came up with an idea: we can set 
up some typical cases as the starting point for our system. For example, if we 
roughly know what kind of problems that the patient may have, for instance, 
tricuspid stenosis, then we can set a tricuspid case as our original parameters 
instead of always starting from normal values. It will certainly save lots of 
iterations and decrease the possibility of getting wrong answers. An intuitive 
explanation to this idea is that again, the distance from the “target” to the 


“original” state is closer. 


e Since clinically ECG is very easy to record and get, we can connect our sys- 
tem with ECG analysis that may help our search process. For example, ECG 
diagnoses for myocardial infarction (MI), hypertrophy, etc., can be used in es- 


tablishing the initial diagnosis. 


e We can also try other ways to do the reasoning process. For example, in our 
current system, we use a flow chart to describe the reasoning process. It is 
part of the reason that most of the answers seem to favor left ventricular prob- 
lem. A change that we can make is that for each iteration, we find the biggest 
difference among all observable data and think about how to decrease this dif- 
ference instead of always starting from matching cardiac output, pulmonary 
wedge pressure, etc. Or we can impose a probabilistic value to each decision 
branch. In that way, we may increase the possibility of escaping from a local 


minimum state. 


e We noticed that part of the reason why sometimes the system seems to get 
stuck, is because it is hard for the system to figure out that, “Hey, I have 
been here before”. How to detect and describe this situation? We may define 
several domains in the multi-dimensional space to categorize the states of the 
system. Once we find out the system is back to a previous state, we can apply 
another reward strategy to increase the unit step size in this situation, so that 


the system can possibly get out from the local minimum state. 
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e In the long term, we hope to apply these methods to actual patient data. More 


summative evaluation methods will then need to be devised. 


e Another problem we noticed is that the algorithm always tried to minimize 
the output errors for only the next step. That means our program lacks a 
global view of the problem. Just like playing the Chess game, sometimes, we 
may have to take a sacrifice in order to win the game eventually. Currently, our 
algorithm is not smart enough to sacrifice the output errors in certain iterations 
to eventually decrease the final total output error. How to make our system 
smarter? We may need to apply some general guidelines for each case, so that 
the system may seem to know which parameter(s) will need to be focused on 


globally. 


6.3 Summary 


In this project, we have developed an iterative method for quantitatively or semi- 
quantitatively analyzing ICU hemodynamic data. The key idea embodied in our 
method is that with the help of a model, it is possible to solve a multi-dimensional 
search problem by iteratively running the model, that is, test and see. The other 
qualitative description may help the search process. As concluded from the above 
discussion, we still have a long way to go before we can fully, accurately and efficiently 
automate the interpretation process. However, we believe that we have developed a 
successful prototype of the final system that will guide us towards our objective. Just 
as the IBM Deep Blue computer has won the Chess game, we are also optimistic that 


our system can be improved for practical use in ICU areas in the foreseeable future. 
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Appendix A 


Case Study Results for Testing 
Expert System 


In this appendix, we will present the case study results for testing our expert syetm. 


For each case, two tables: 


e Comparison of Pseudo-Patient Input Parameters with Estimated Input Param- 


eters 
e Comparison of Output Variables of Pseudo-Patient with those of Simulator 


and two figures: 


e Total Output Error Function vs. Number of Iterations 


e Total Input Error Function vs. Number of Iterations 


will be presented. 
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A.1 Cases Without Compensatory Effects 


A.1.1 Blood Loss 


Psendo-Patient 


Heart Rate (beats/min) 


Total Blood Volume (ml) 5000 [| 4500's 
Trans-thoracic Pressure (nmHg) 


Capacitances: (ml/mmHg) 


Table A.1: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Blood Loss Case 


Measurable Variables 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) 
Cardiac Output (ml/min) 

Heart Rate (beats/min) 

Mean Right Ventricular Pressure (mmHg) 
Mean Pulmonary Arterial Pressure (mmHg) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 
Total Number of Iterations 


Table A.2: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Blood Loss Case 


lThe items in bold type are different from the normal values. 
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Blood Loss — 


3,5 


2.5 


1,5 


0,5 
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1 2 3 4 5 6 7 


Figure A-1: Input Error Function vs. Number of Iterations in Blood Loss Case 


Blood Loss —— 
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Figure A-2: Output Error Function vs. Number of Iterations in Blood Loss Case 


oo 


oO 


A.1.2 Arteriolar Dilation 


Pseudo-Patient 


Heart Rate (beats/min) 


Total Blood Volume (ml) 5000 5000 
Trans-thoracic Pressure (mmHg) 


Capacitances: (ml/mmHg) 


Table A.3: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Arteriolar Dilation Case 


Measurable Variables 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) 
Cardiac Output (ml/min) 

Heart Rate (beats/min) 

Mean Right Ventricular Pressure (mmHg) 
Mean Pulmonary Arterial Pressure (mmHg) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 


Total Number of Iterations Total Output Error 


Table A.4: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Arteriolar Dilation Case 
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Arteriolar Dilation — 


0 0,5 1 1.5 2 2.5 3 3,5 4 
Figure A-3: Input Error Function vs. Number of Iterations in 
Arteriolar Dilation Case 


Arteriolar Dilation — 


0 0,5 1 1.5 2 2,5 3 3,5 4 
Figure A-4: Output Error Function vs. Number of Iterations in 
Arteriolar Dilation Case 
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A.1.3 Left Ventricular Systolic Dysfunction 


Psendo-Patient 


Heart Rate (beats/min) 


Total Blood Volume (ml) 5000 5000 
Trans-thoracic Pressure (mmHg) 


Capacitances: (ml/mmHg) 


Table A.5: Comparison of Pseudo-Patient Input Parameters with Estimated 
Input Parameters in Left Ventricular Systolic Dysfunction Case 


Measurable Variables 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) 
Cardiac Output (ml/min) 

Heart Rate (beats/min) 

Mean Right Ventricular Pressure (mmHg) 
Mean Pulmonary Arterial Pressure (mmHg) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 


Total Number of Iterations Total Output Error 


Table A.6: Comparison of Output Variables of Pseudo-Patient with those of 
Simulator in Left Ventricular Systolic Dysfunction Case 
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Left Ventricular Systolic Dysfunction — 


Figure A-5: Input Error Function vs. Number of Iterations in 
Left Ventricular Systolic Dysfunction Case 


Left Ventricular Systolic Dysfunction — 


Figure A-6: Output Error Function vs. Number of Iterations in 
Left Ventricular Systolic Dysfunction Case 
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A.1.4 Arteriolar and Venous Dilation 


Pseudo-Patient 


Heart Rate (beats/min) 


Total Blood Volume (ml) 5000 5000 
Trans-thoracic Pressure (mmHg) 


Capacitances: (ml/mmHg) 


ce 


Table A.7: Comparison of Pseudo-Patient Input Parameters with Estimated 
Input Parameters in Arteriolar and Venous Dilation Case 


Measurable Variables 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) 
Cardiac Output (ml/min) 

Heart Rate (beats/min) 

Mean Right Ventricular Pressure (mmHg) 
Mean Pulmonary Arterial Pressure (mmHg) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 


Total Number of Iterations | 9 | Total Output Error 


Table A.8: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Arteriolar and Venous Dilation Case 
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Arteriolar and Venous Dilation — 


0 1 2 3 4 5 6 7 8 
Figure A-7: Input Error Function vs. Number of Iterations in 
Arteriolar and Venous Dilation Case 


Arteriolar and Venous Dilation — 


Figure A-8: Output Error Function vs. Number of Iterations in 
Arteriolar and Venous Dilation Case 
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A.2 Cases With Compensatory Effects 


A.2.1 Hypertension 


Teart Rate (beats/min) 
[Total Blood Volume (mly | 5000 | _5000_| 5268.55 |_| 
Capacitances: (ml/mmHg) 
Pp WaystSSSSCSC~dYSCi CdCI 
Wiest | 100 [6.0 | 6.07 | 1.2%_| 
[Arenal | 16 [0.9 [_0.a6 | 4.4%_| 
Zero-Pressure Volumes: (ml) 
[Venous 2500.0 [2200.0 | 2500.0 | 
[Pais are 90.0 90.0980. 
Resistances: (mmHg*sec/ml) 
[Microvascular | 1.0 | 2.00 | 1.86 [2% | 


Table A.9: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Hypertension Case 


Mean Central Venous Pressure (mmHg) | _66_| 74 | __73 | 14% | 


Table A.10: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Hypertension Case 
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Hypertension —— 


3,5 


2.5 


1.5 


0,5 


0 
0 5 10 15 20 25 


Figure A-9: Input Error Function vs. Number of Iterations in Hypertension Case 


Hypertension — 


3,5 


2.5 


1,5 


0,5 


0 5 10 15 20 25 


Figure A-10: Output Error Function vs. Number of Iterations in Hypertension Case 
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A.2.2 Right Ventricular Infarction 


[Total Blood Volume (mly_| 5000] 5000 | 5000] 
Capacitances: (ml/mmHg) 

0 
PRvayst SSS «dCi 
[RV diast_—— 0.0 8 2 | 


4 
.40 


Pum, Ard 910 90.0 90.0 | 0% | 
Resistances: (mmHg*sec/ml) 

[Pulmonary | 0.08_[-_ 0.08 | 0.09 | 12.5% 
[Total Number ofTrerations | oc? | Total Input Error SBM 


Table A.11: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Right Ventricular Infarction Case 


Error 


55 
62 


Total Number of Iterations Total Output Error 


Table A.12: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Right Ventricular Infarction Case 


250 means the program did not stop by itself. 
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Right Ventricular Infarction — 


0 50 100 150 200 250 
Figure A-11: Input Error Function vs. Number of Iterations 
in Right Ventricular Infarction Case 


Right Ventricular Infarction — 


0 50 100 150 200 250 
Figure A-12: Output Error Function vs. Number of Iterations 
in Right Ventricular Infarction Case 
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A.2.3. Hypertrophic Cardiomyopathy 


Psendo-Patient 


Heart Rate (beats/min) 
Total Blood Volume (ml) 


[10.0 [3.00 
172 
[00 [_4.0__—i ata 


15.0 15.0 


490.0 490.0 490.0 
1175 2675 2740.28 


Table A.13: Comparison of Pseudo-Patient Input Parameters with Estimated 
Input Parameters in Hypertrophic Cardiomyopathy Case 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) | 66 [| 204 [ 206 | 
Cardiac Outpat (ml/min) 
Heart Rate (beats/min) 100 
Mean Right Ventricular Pressure (mmHg) 


Mean Pulmonary Arterial Pressure (mmHg) 
Maximum Pulmonary Arterial Pressure (mmHg) 
Minimum Pulmonary Arterial Pressure (mmHg) 
Pulmonary Wedge Pressure (mm) 
Teft Heart Ejection Fraction (72) 
Hight Heart Fjection Fraction (%) 


Table A.14: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Hypertrophic Cardiomyopathy Case 
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Hypertrophic Cardiomyopathy — 


0 10 20 30 40 50 60 70 
Figure A-13: Input Error Function vs. Number of Iterations 
in Hypertrophic Cardiomyopathy Case 


Hypertrophic Cardiomyopathy —— 


0 10 20 30 40 50 60 70 
Figure A-14: Output Error Function vs. Number of Iterations 
in Hypertrophic Cardiomyopathy Case 
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80 


80 


A.2.4 Pulmonary Hypertension 


Psendo-Patient 


Heart Rate (beats/min) 
Total Blood Volume (ml) 


Table A.15: Comparison in Pseudo-Patient Input Parameters with Estimated 
Input Parameters in Pulmonary Hypertension Case 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) | 66 [| 15 [120 | 
Cardiac Outpat (ml/min) 
Heart Rate (beats/min) 120 120 

Mean Right Ventricular Pressure (mmHg) 


Mean Pulmonary Arterial Pressure (mmHg) 
Maximum Pulmonary Arterial Pressure (mmHg) 
Minimum Pulmonary Arterial Pressure (mmHg) 
Pulmonary Wedge Pressure (mam) 
Tett Heart Ejection Fraction (72) 
Hight Heart Fjeotion Fraction (%) 
[Total Number of erations | oo | Total Output Error 


Table A.16: Comparison of Output Variables of Pseudo-Patient 
with those of Simulator in Pulmonary Hypertension Case 


gg 


Pulmonary Hypertension —— 


0 50 100 150 200 250 300 
Figure A-15: Input Error Function vs. Number of Iterations 
in Pulmonary Hypertension Case 


Pulmonary Hypertension —— 


0 50 100 150 200 250 300 
Figure A-16: Output Error Function vs. Number of Iterations 
in Pulmonary Hypertension Case 


100 


A.2.5 Large AV Fistula 


Psendo-Patient 


Heart Rate (beats/min) 
Total Blood Volume (ml) 


[0.08 [0.04 


Table A.17: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in AV Fistula Case 


Measurable Variables 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) 
Cardiac Output (ml/min) 

Heart Rate (beats/min) 

Mean Right Ventricular Pressure (mmHg) 
Mean Pulmonary Arterial Pressure (mmHg) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 


Total Number of Iterations Total Output Error 


Table A.18: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in AV Fistula Case 
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AV Fistula — 
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Figure A-17: Input Error Function vs. Number of Iterations in AV Fistula Case 


AV Fistula — 


3,5 


2.5 


1.5 


0,5 
0 


0 5 10 15 20 25 


Figure A-18: Output Error Function vs. Number of Iterations in AV Fistula Case 
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A.2.6 Dilated Cardiomyopathy and Bi-Ventricular Conges- 


tive Heart Failure 


Pseudo-Patient 


Heart Rate (beats/min) 
Total Blood Volume (ml) 


11.0 


430d 


8.40 


0. 
[2500.0 [2200 


[10 [0.9 if 0.88 


Table A.19: Comparison of Pseudo-Patient Input Parameters with Estimated 
Input Parameters in Dilated Cardiomyopathy and 
Bi-Ventricular Congestive Heart Failure Case 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) 
Cardiac Output (ml/min) 


Heart Rate (beats/min) 125 
Mean Right Ventricular Pressare (mmlTg) 


Mean Pulmonary Arterial Pressure (mmHg) 
Maximum Pulmonary Arterial Pressure (mmHg) 
Minimum Pulmonary Arterial Pressure (mmHg) 
Pulmonary Wedge Pressure (mmllg) 
Teft Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 
Total Number of Terations 


Table A.20: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Dilated Cardiomyopathy and 
Bi-Ventricular Congestive Heart Failure Case 
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Dilated Cardiomyopathy with Bi-Ventricular Congestive Heart Failure —— 


0 10 20 30 40 50 60 70 80 

Figure A-19: Input Error Function vs. Number of Iterations in 
Dilated Cardiomyopathy and Bi-ventricular 
Congestive Heart Failure Case 


Dilated Cardiomyopathy with Bi-Ventricular Congestive Heart Failure —— 


30 


0 10 20 30 40 50 60 70 80 

Figure A-20: Output Error Function vs. Number of Iterations in 
Dilated Cardiomyopathy and Bi-ventricular 
Congestive Heart Failure Case 
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30 


A.2.7 Acute Myocardial Infarction with Bi- Ventricular Heart 


Failure 


Pseudo-Patient 


Heart Rate (beats/min) 
Total Blood Volume (ml) 


[30 [2.0 


8.40 


Table A.21: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Acute Myocardial 
Infarction and Bi-ventricular Heart Failure Case 


Measurable Variables 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) 
Cardiac Output (ml/min) 

Heart Rate (beats/min) 

Mean Right Ventricular Pressure (mmHg) 
Mean Pulmonary Arterial Pressure (mmHg) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 


Table A.22: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Acute Myocardial Infarction 
and Bi-Ventricular Heart Failure Case 
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Acute Myocardial Infarction with Bi-Ventricular Heart Failure —— 


0 10 20 30 40 50 60 70 80 
Figure A-21: Input Error Function vs. Number of Iterations in 
Acute Myocardial Infarction and Bi-ventricular Heart Failure Case 


Acute Myocardial Infarction with Bi-Ventricular Heart Failure —— 


0 10 20 30 40 50 60 70 80 
Figure A-22: Output Error Function vs. Number of Iterations in 
Acute Myocardial Infarction and Bi-Ventricular Heart Failure Case 
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A.2.8 Septic Shock 


Pseido-Patient 


Heart Rate (beats/min) 
Total Blood Volume (ml) 


Table A.23: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Septic Shock Case 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) | 66 [| 39 [ 40 | 
Cardiac Outpat (ml/min) 
Heart Rate (beats/min) 150 150 

Mean Right Ventricular Pressure (mmHg) 


Mean Pulmonary Arterial Pressure (mmHg) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 


Table A.24: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Septic Shock Case 
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Septic Shack —— 
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Figure A-23: Input Error Function vs. Number of Iterations in Septic Shock Case 


Septic Shack —— 


3,5 


2.5 


1,5 


0,5 


0 10 20 30 40 50 60 


Figure A-24: Output Error Function vs. Number of Iterations in Septic Shock Case 
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A.2.9 Pulmonary Hypertension with Right Heart Failure 


Heart Rate (beats/min) 
Total Blood Volume (ml) } 5000 [| 6500 | 
Trans-thoracic Pressure (mmHg) 


Capacitances: (ml/mmHg) 


Table A.25: Comparison of Pseudo-Patient Input Parameters with Estimated Input 
Parameters in Pulmonary Hypertension and Right Heart Failure Case 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) | 66 [| 228 [221 | 
Cardiac Outpat (ml/min) 
Heart Rate (beats/min) 
Mean Right Ventricular Pressure (mmHg) 


Mean Pulmonary Arterial Pressure (mmHg) 
Maximum Pulmonary Arterial Pressure (mmHg) 
Minimum Pulmonary Arterial Pressure (mmHg) 
Pulmonary Wedge Pressure (mm) 
Teft Heart Ejection Fraction (72) 
Hight Heart Ejection Fraction (7) 


Table A.26: Comparison of Output Variables of Pseudo-Patient with those of 
Simulator in Pulmonary Hypertension and Right Heart Failure Case 
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Pulmonary Hypertension with Right Ventricular Heart Failure — 


0 5 10 15 20 25 30 35 40 45 50 
Figure A-25: Input Error Function vs. Number of Iterations in 
Pulmonary Hypertension and Right Heart Failure Case 


Plilonary Hypertension with Right Ventricular Heart Failure — 


0 5 10 15 20 25 30 35 40 45 50 
Figure A-26: Output Error Function vs. Number of Iterations in 
Pulmonary Hypertension and Right Heart Failure Case 


110 


A.2.10 Mitral Stenosis I 


Pseudo-Patient 


Heart Rate (beats/min) 
Total Blood Volume (ml) 


15.0 


490.0 490.0 490.0 
1175 1675 1705.11 


A 
[0.08 [0.80] 


Table A.27: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Mitral Stenosis Case I 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) | 66 [| 78 [ 82 | 
Cardiac Output (ml/min) 
Heart Rate (beats/min) 100 


Mean Right Ventricular Pressure (mig) 
Mean Pulmonary Arterial Pressure (mmFig) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 


Table A.28: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Mitral Stenosis Case I 


Lact 


Mitral Stenosis I — 


0 5 10 15 20 25 
Figure A-27: Input Error Function vs. Number of Iterations 
in Mitral Stenosis Case I 


Mitral Stenosis I — 


0 5 10 15 20 25 
Figure A-28: Output Error Function vs. Number of Iterations 
in Mitral Stenosis Case I 


1d? 


30 


30 


A.2.11 Mitral Stenosis II 


Teart Rate (beats/min) 
[Total Blood Volume (ml) | 5000 | _5500___| 844.85 _[ | 
Capacitances: (ml/mmHg) 
Pp Layst SS SCSCSC~dYSC CSCC 
T2.5% 
P_RVayst_——SCSCSCSCS 20 i dCi |_| 
Zero-Pressure Volumes: (ml) 
[Pam are 90.0 90.0 90.0 J 
Resistances: (mmHg*sec/ml) 
[Winfow—SSSCSC~iSCi | Cid SC 
[Pulmonary | 0.08 | ~_0.80 | 0.29 [3.3%] 


Table A.29: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Mitral Stenosis Case II 


Mean Central Venous Pressure (mmHg) | _66_| 78 | __7®8 | 0% 


Table A.30: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Mitral Stenosis Case II 
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Mitral Stenosis I] — 


0 5 10 15 20 25 30 
Figure A-29: Input Error Function vs. Number of Iterations 
in Mitral Stenosis Case II 


Mitral Stenosis I] — 


0 5 10 15 20 25 30 
Figure A-30: Output Error Function hvjkcxvs. Number of Iterations 
in Mitral Stenosis Case II 
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A.2.12 Tricuspid Stenosis 


Heart Rate (beats/min) 
Total Blood Volume (ml) } 5000 [| 5500 | 
Trans-thoracic Pressure (mmHg) 


Capacitances: (ml/mmHg) 


Table A.31: Comparison of Pseudo-Patient Input Parameters with 
Estimated Input Parameters in Tricuspid Stenosis Case 


Mean Systemic Arterial Pressure (mmHg) 
Maximum Systemic Arterial Pressure (mmHg) [ 113.2 [ 955 [96.0 
Minimum Systemic Arterial Pressure (mmHg) 
Mean Central Venous Pressure (mmHg) | 66 [| 123.0 [132 | 
Cardiac Output (ml/min) 
Heart Rate (beats/min) 120 
Mean Right Ventricular Pressure (mmHg) 


Mean Pulmonary Arterial Pressure (mmMg) 


Pulmonary Wedge Pressure (mmHg) 
Left Heart Ejection Fraction (%) 
Right Heart Ejection Fraction (%) 


Table A.32: Comparison of Output Variables of Pseudo-Patient with 
those of Simulator in Tricuspid Stenosis Case 
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Tricuspid Stenosis —— 


0 5 10 15 20 
Figure A-31: Input Error Function vs. Number of Iterations 
in Tricuspid Stenosis Case 


Tricuspid Stenosis —— 


0 5 10 15 20 
Figure A-32: Output Error Function vs. Number of Iterations 
in Tricuspid Stenosis Case 
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